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P\J , Abstract Starting from one-range addition theorems for Slater-type functions, which 

are expansion in terms of complete and orthonormal functions based on the general- 
ized Laguerre polynomials, Guseinov constructed addition theorems that are expan- 
sions in terms of Slater-type functions with a common scaling parameter and inte- 
gral principal quantum numbers. This was accomplished by expressing the complete 
and orthonormal Laguerre-type functions as finite linear combinations of Slater-type 
functions and by rearranging the order of the nested summations. Essentially, this 
1-^ ■ corresponds to the transformation of a Laguerre expansion, which in general only 

converges in the mean, to a power series, which converges pointwise. Such a trans- 
it ■ formation is not necessarily legitimate, and this contribution discusses in detail the 
jrt I difference between truncated expansions and the infinite series that result in the ab- 
sence of truncation. 
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1 Introduction 

Electronic structure theory is a highly interdisciplinary research topic which benefited 
greatly from interactions with related scientific disciplines. It is generally agreed that 
molecular electronic structure calculations only became feasible because of the spec- 
tacular advances in computer hard and software. But we must not forget that advances 
in pure and applied mathematics also played a crucial role. 

While mathematicians and theoretical physicists have a venerable tradition of a 
close and mutually beneficial collaboration, the same cannot be said about the inter- 
action of mathematicians and theoretical chemists (my personal views on this topic 
are explained in [208]). This is highly deplorable. Theoretical chemists could - and 
should - learn more about new mathematical concepts or powerful numerical tech- 
niques. On the other hand, electronic structure theory offers quite a few challenging 
mathematical and computational problems which could serve as a valuable source of 
inspiration for mathematicians. A more extensive collaboration would certainly be 
mutually beneficial (see for example ||208[ Section 8] and references therein). 

Advances in mathematics are particularly important for molecular multicenter 
integrals, which occur in molecular calculations on the basis of the Hartree-Fock- 
Roothaan equations JSOl Il36l Il67ll and also in other approximation schemes. The 
efficient and reliable evaluation of the three- and six-dimensional molecular inte- 
grals, which are notoriously difficult in the case of the physically better motivated 
exponentially decaying basis functions, is one of the oldest mathematical and com- 
putational problems of molecular electronic structure theory. In spite of heroic efforts 
of numerous researchers, no completely satisfactory solution has been found yet, and 
we have to continue searching for more powerful mathematical and computational 



techniques. A review of the older literature can be foun d in a rticles by Browne 1118 1 
Dalgarno 0230 . Harris and Michels 1 135], and Huzinaga 114011 . 



Quantum mechanics only determines which types of molecular integrals we have 
to evaluate, but how we manipulate and ultimately evaluate them is a mathematical 
problem. Thus, for a researcher working on multicenter integrals, physical insight and 
a good knowledge of quantum mechanics is actually less important than the ability 
of skillfully manipulating complicated expressions involving special functions and a 
profound knowledge of advanced mathematical techniques with a special emphasis 
on numerics. 

The evaluation of multicenter integrals is difficult because of the different centers 
occurring in their integrands. This effectively prevents the straightforward separa- 
tion of the three- and six-dimensional integrals into products of simpler integrals. 
A promising computational strategy requires that we find a way of separating the 
integration variables at tolerable computational costs. 

Principal tools, which can accomplish a separation of integration variables, are so- 
called addition theorems. These are special series expansions of a function f{r±r') 
with r,r' E M? in terms of other functions that only depend on either r or r'. The 
basic features of these fairly complicated series expansions are reviewed in Section 
121 A much more detailed treatment will be given in my forthcoming review |211]. 

Infinite series expansions are among the most fundamental mathematical tools 
with countless applications not only in mathematics, but also in science and engi- 



neering. Nevertheless, a mathematically rigorous use of infinite series is not neces- 
sarily an easy thing. Many scientists have a highly cavalier attitude when it comes 
to questions of existence and convergence. Since mathematics is ultimately used to 
describe natural phenomena, it is tempting to believe that nature guarantees that all 
intermediate mathematical manipulations are legitimate. Needless to say that such an 
attitude is overly optimistic and can easily have catastrophic consequences. 

Strictly speaking, an infinite series is a meaningless object unless we specify its 
convergence type - for example pointwise convergence, convergence in the mean, 
or even distributional or weak convergence - and provide convincing evidence that 
this series converges according to its specified type of convergence. It is important 
to take into account that different convergence types of series expansions imply dif- 
ferent mathematical properties with different advantages and disadvantages. This ap- 
plies also to addition theorems. As discussed in my forthcoming review 121 ill , the 
different convergence types of addition theorems may well provide the most useful 
characterization of their properties. 

Slater-type functions were originally introduced by Slater llTSl 117911 to provide 



computationally convenient analytical approximations to numerically determined so- 
lutions of effective one-particle Schrodinger equations. In unnormalized form, they 
can be expressed as follows: 

X^,Lil^,'-) = {firf-^-'^-P'-^J^{M = (^rf-'e-/^'-yf (r/r). (1.1) 

Here, j3 > is a scaling parameter, '3^1' {fir) — {j5r)^Y^ (r/r) is a regular soHd hai"- 
monic and Y^{r/r) is a surface spherical harmonic (in my own work, I have always 
used the phase condition of Condon and Shortley [21, Eqs. (6) and (9) on p. 115]). 
By a slight abuse of language, the index A^ is frequently called principal quantum 
number In most, but not in all cases A^ is a positive integer satisfying N >L+\. 

Because of the importance of Slater-type functions as basis functions in atomic 
and molecular electronic structure calculations, it is not surprising that there is an 
extensive literature on their notoriously difficult multicenter integrals in general as 
well as on their addition theorems in special. A reasonably complete bibliography on 
their addition theorems would be beyond the scope of this article. Let m e just mention 



two classic articles by Barnett and Coulson [^ and by Lowdin 1115311 . respectively, 
which have inspired many other researchers. 



As discussed in more detail in Section[3j Guseinov 14 ll - l44, 14611 derived so-called 
one-range addition theorems for Slater-type functions with integral and nonintegral 
principal quantum numbers. He expanded Slater-type functions XfjiiPy''^''') ^^ 
terms of the functions { k^"\{7^ '') } ^ ' which are complete and orthonormal in cer- 
tain Hilbert spaces and which are defined in ( 12.131 ). The radial parts of these functions 
are based on the generalized Laguerre polynomials whose most inportant properties 
are reviewed in AppendixiDl 

On the basis of the approach developed in f4 1-44, 461, Guseinov and coworkers 
produced an amazing number of articles on addition theorems and related topics 1401 
l45[l47U49ll5lU53ll5dJ74l76 Ul03lll07[ll08[llllUll9lll2lU l34^. 

The immediate reason for writing this article is a one-center expansion con- 



structed and applied by Guseinov and Mamedov 1112211 . They expanded a Slater- 
type function XNiiP^'') with an in general nonintegral principal quantum number 



A^ e M \ N in terms of Slater-type functions {z^L(7j'")}°°-f + 1 with integral princi- 
pal quantum numbers « G N and an in general different common scaling parameter 
77^ j3 > 1 122, Eq. (4)]. This expansion was used by Guseinov and Mamedov for the 
construction of what they called series expansions for overlap integrals of Slater-type 
functions with nonintegral principal quantum numbers in terms of overlap integrals 
of Slater-type functions with integral principal quantum numbers. 

For Guseinov and coworkers, who had done a lot of work on Slate r-typ e functions 
with nonintegral principal quantum number s 1451 ISOl 52l 54 55l I6II |68l |72|, |75l [77 , 
io, 82-84, 90, 91, 93-96, 98, 102-106, 109- (lHfriin23ni3lll . such an expansion 



would be extremely useful: Slater-type functions with nonintegral principal quantum 
numbers could be replaced by Slater-type functions with integral principal quantum 
numbers whose multicenter integrals can be evaluated (much) more easily. 

It is, however, trivially simple to show that such a one-center expansion, which in 
my notation can be written as follows, 

X^,L{P,r)= f Xr^'(j3,7)Z„'l(7,'-), N eR, (1.2) 

n=L+l 

only exists if the principal quantum number A^ is a positive integer satisfying A^ > 
L+1. If A^ is nonintegral or zero, i.e., if A^ G K \ N, this expansion does not exist. 



The nonexistence of this expansion also played a major role in my discussion 11207 



Section 7] of Guseinov's treatment of one-range addition theorems for Slater-type 



functions. However, my discussion in 112071 Section 7] was incomplete. Its scope is 
extended considerably by this article. 

From a methodological point of view, we only have to utilize the obvious fact 
that an expansion of XmlH^t'') ^^ terms of Slater-type functions {x^LiYi'')} -r+i 
with integral principal quantum numbers n and a common scaling parameter 7 > 
is nothing but a power series expansion of exp(7r);if^^()3 , r) about r = in disguise. 
This interpretation of the expansion (11.21 ) offers some valuable insight. Relatively 
little is known about the convergence and existence of expansions in terms of Slater- 
type functions {z,fL(7, '")} , but it is normally comparatively easy to decide whether 
such a corresponding power series expansion about r = exists and for which values 
of r it converges. 

We can convert the expansion (ll.2l i to a more transparent one-dimensional radial 
problem by canceling the spherical harmonics. Then, the right-hand side of ( ll.2l i can 
be converted to a power series in y by multiplying either side of ( 11.2b by exp(7r); 

e(r-/3y(^,)iv-i = £ Xf^)(^,7)(7'-r'. (1.3) 

n=L+l 

Let us now assume that the principal quantum number A^ is a positive integer 
satisfying A^ > L + 1 . Then we only have to replace the exponential on the left-hand 
side of (11.3b by its power series to obtain: 

V=0 ^- n=L+l 



Comparison of the coefficients of equal powers yields explicit expressions for the 
coefficients X„ ' (jS,/). 

But if the principal quantum number A^ is not a positive integer, the matching of 
the coefficients of equal powers does not work: On the left-hand side of (ll.4l l there 
are either only nonintegral powers or some negative powers, and on the right-hand 
side there are only integral and nonnegative powers. Accordingly, a power series 
expansion for e^''^^''' (/3r)^^' can only exist if A^ is a positive integer satisfying A^ > 
L+1. 

Sin ce ev ery power series is also a Taylor series for some function (see for ex- 
ample il55h . the expansion (11.2b can also be derived by doing a Taylor expansion 
of exp((7~ j3)r))(/3r)^^' about r — 0, yielding the same conclusions. Fractional 
or nonintegral powers r" with a G M \ No and thus also exp([7— j3]r)r'^^' with 
A' e K \ N do not have continuous derivatives of arbitrary order at r = 0. Thus, 
A^ e K \ N implies that the leading expansion coefficients X), ' (/3 , 7) with indices 
n < [N\ are zero, but all the remaining coefficients with indices n > [N\ are infinite 
in magnitude. Here, [A^J stands for the integral part of A^. 

The nonexistence problems mentioned above are a direct consequence of the ob- 
vious fact that the radial part of a Slater-type function with a nonintegral or negative 
principal quantum number is not analytic in the sense of complex analysis at r = 0. 
This is so elementary that it is hard to understand why nobody had noticed the nonex- 
istence of the expansion (ll.2l i for A^ G R \ N before. 

In spite of the nonexistence of (11.2b for A^ G M \ N, Guseinov and Mamedov 



012211 had used their version of this one-center expansion for numerical purposes 
and presented apparently meaningful numerical results for overlap and other, closely 
related integrals of Slater-type functions with nonintegral principal quantum numbers 



in III22I Tables 1 and 2] . 

If these numerical results are genuine - which I assume - then there is only one 
log icall y satisfactory conclusion: Contrary to their claim (see the text before Eq. (4) 



of 1112211 '). Guseinov and Mamedov did not use expansions involving an infinite num- 
ber of Slater-type functions with integral principal quantum numbers, which - as 
discussed above - do not exist for A^ G M \ N. Instead, they only employed some 
approximations consisting of ?i finite number of terms and extrapolated - albeit incor- 
rectly - that these finite approximations remain meaningful in the limit of an infinite 
number of terms. 

As is well known from the mathematical literature, truncated expansions, which 
are approximations consisting of a finite number of terms, may have mathematical 
properties that are not at all related to those of a complete expansion consisting of 
an infinite number of terms. In particular, it can happen that a truncated expansion 
can be mathematically meaningful as well as numerically useful, although its limit 
as an expansion of infinite length does not exist (see also Appendix IE] where the 
semiconvergence of certain infinite series is reviewed briefly). It is the purpose of 
this article to investigate this as well as some closely related question. 

The one-center expansion used by Guseinov and Mamedov ['l22, Eq. (4)] was 
originally derived by Guseinov [46, Eq. (21)] as the one-center limit of a class of 
addition theorems for Slater-type functions [46, Eq. (15)]. Guseinov's approach is 



based on one-range addition theorems for Slater-type functions, which are expan- 
sions in terms of Guseinov's Laguerre-type functions {k^"\{'Y, '') } c ■ These func- 
tions are then replaced by Slater-type functions according to ( 13.21 ). Ultimately, this 
yields addition theorems that are expansions in terms of Slater-type functions with 
integral principal quantum numbers. Unfortunately, the derivation of these expan- 
sions is fa i rly d ifficult and involves manipulations whose validity is not at all obvious 



l 204lE()6ll207ll . 



In view of the amazing number of articles on addition theorems and related topics 
pubUshedby Guseinov and coworkers |40-49, 51-53, 56-74, 76-103, 107, 108, lllj- 
119lll2ll4l34ll . such an investigation should be of some principal interest. 



This article sheds some light on certain more subtle differences of orthogonal 
and nonorthogonal expansions which - although in principle known - are often not 
sufficiently appreciated and taken into account in the literature on electronic structure 
calculations. 

In Section 121 basic features of both one-range and two-range addition theorems 
are reviewed. Section [3] discusses Guseinov's derivation of expansions in terms of 
Slater-type functions by rearranging his one-range addition theorems for Slater-type 
functions with a special emphasis on numerical stability issues. From a mathematical 
point of view, Guseinov's treatment of addition theorems ultimately boils down to 
the transformation of a Laguerre series to a power series. The most impor tant p rop- 
erties of this transformation, whose theory had recently been formulated in 120711 . are 
reviewed in Section|4l 

The central results of this manuscript are presented in Sections|5]-|2] In Section|5] 
it is shown that the transformation formulas described in Section |4] fully explain the 
nonexistence of Guseinov's rearrangements of one-center expansions in the case of 
nonintegral principal quantum numbers. 

In the case of the simpler one-center expansions, which correspond to the La- 
guerr e series ( 15.21 ) for z^ with p G M \ Nq, this analysis had already been done in 



1 2071 Section 3]. But my results for the more complicated one-center expansions. 



whic h correspond to the Laguerre series ( 15.11 ) for z^e""' with p e M \ No, are new. In 



1 2071 Section 7] it was only stated on the basis of general analyticity principles that 
the Laguerre series (15.11 1 cannot exist if p G M \ No, but a detailed analysis of the 
transformation formulas could not be done. At that time I did not know yet the large 
index asymptotics of some special Gaussian hypergeometric series 2F\ which were 
derived in Section |5] These asymptotic expressions make an asymptotic analysis of 
the transformation formulas possible. 

As discussed in more details in Section |2l a one-range addition theorem for a 
function /(r ± r') provides for all arguments r, r' G M^ a unique representation of / 
in separated form that is valid for the whole argument set M? x M?. Since the one- 
center limits of Guseinov's rearranged one-range addition theorems do not exist in 
the important special case of nonintegral principal quantum numbers, it was argued 



in 112071 Section 7] that at least some of Guseinov's rearrangements cannot be one- 
range addition theorems. How ever, the exact nature of Guseinov's rearrangements 
was not yet understood in l'207'l. 

Because of technical problems, a straightforward analysis of the transformation 
formulas described in Section |4] is not possible in the case of Guseinov's one-range 



addition theorems. However, in Section |6] it is shown that the rearrangement of Gu- 
seinov's one-range addition theorems produces two-range addition theorems. This 
follows at once from the fact that a Slater-type function XnlH^^''^ '"') '^^ equiva- 
lently the function s^' Xnl iP'''^''') i^ singular for r ± r' = 0. 

Section|7]discusses some numerical implications of the results derived in Sections 
|5]and|6l The most important result is that rearrangements of Guseinov's truncated 
one-center expansions are semiconvergent with respect to a variation of the truncation 
order if the principal quantum number is nonintegral. This follows from the leading 
order asymptotic approximation ( 17.21 ). 

This article is concluded by a summary and an outlook in Section[8] Appendices 
lAl - lEl review general aspects of series expansions, the largely complementary fea- 
tures of power series and orthogonal expansions, the for our purposes most important 
properties of generalized Laguerre polynomials, and semiconvergent expansions. 



2 Basic Features of Addition Theorems 

Addition theorems are series expansions of a function f{r ± r') with r, r' G K.^ in 
terms of other functions that depend only on either r or r'. Because of the complete- 
ness and the convenient orthogonality properties of the (surface) spherical harmonics 
Y"\9,^), it is customary to express the angular parts of addition theorems in terms 
of spherical harmonics whose arguments correspond to the solid angles r/r and r'/r', 
respectively. 

In atomic and molecular electronic structure theory, we are predominantly inter- 
ested in addition theorems for functions 

Frir)=Mr)Yrir/r), (2.1) 

which can be factored into a radial part f({r) multiplied by a surface spherical har- 
monic. In the language of angular momentum theory, such a function F'"{r) is an 
irreducible spherical tensor with a fixed rank £ £ Nq. 

The best known and probably also the most simple addition theorem of such an 
irreducible spherical tensor is the Laplace expansion of the C oulomb or Newton po- 
tential 1 / 1 r ± / 1 . As discussed in Hobson's classic book 111381 §11 on pp. 16-17], the 
mathematical tools needed for the construction of this addition theorem were devel- 
oped by Laplace and Legendre already in the late 18th century. In modern notation, 
the Laplace expansion can be expressed as follows: 



\r±r'\ ■'^21 + 1 r^+i 

|r<| =min(r,r'), |r>| =max(r,r') . (2.2) 

The Laplace expansion is the simplest prototype for a large class of addition the- 
orems that possess a characteristic two-range form. Numerous techniques for the 
derivation of two-range addition theorems are described in the literature. With the 



possible exception of the so-called Fourier transform method advoc ated a lmost si 



multaneously but independently by Ruedenberg 1 168] and Silverstone II177I1 . they are 
all related to power series expansions in one way or the other. 

For example, the Laplace expansion (I2.2l i can be deriv ed ea sily with the help of 
the generating function of the Legendre polynomials P„ [x) \\5A p. 232] : 

[l-2xt + t\^'^ = Y,Pn{x)t" , -1<X<1, |f|<l. (2.3) 

Obviously, this generating function is nothing but a power series expansion about 
t ~Q. For the derivation of the Laplace expansion (12.2b, we only need the so-called 
spherical harmonic addition theorem (see for example il9li Eq. (1.2-3a)]) 

PiicosO) = ^T-y E [}-/'(«/«)] 'J-riv/v), COS0 = — . (2.4) 

The two-range form of the Laplace expansion ( 12.21 ) is necessary to ensure conver- 
gence, because the generating function ( 12.31 ) converges only for |f| < 1. Thus, the 
Laplace expansion converges po/«fwKe for r</r> < 1 and diverges for r</r> > 1. 

The use of generating functions of the Gegenbauer polynomials, which general- 
ize the generating function ( 12.3b of the Legendre polynomials, for the construction of 
addition theorems of spherically symmetric functions will be discussed briefly in Sec- 
tion|6] The two-range forms of these addition theorems follow from the convergence 
conditions of the corresponding generating functions which are nothing but special 
power series expansions, usually in the variable r/r' or more precisely in r^/r^. 

The most principal approach for the construction of pointwise convergent addi- 
tion theorems consists in interpreting such an addition theorem as a three-dimensional 
Taylor expansion (see for example 1 12, p. 181]): 

f{r±r') = £ ^^/(±r') = e-^'/lir') . (2.5) 

H=() 

Thus, the translation operator e*"'^ = ^™l°'^ e-V"' ''■'' t~ I '• generates /(r ± /) by con- 
structing a three-dimensional Taylor expansion of / about ±r' with shift vector r. 
Since the variables r and / are separated, the series expansion (12.5b is indeed an 
addition theorem. 

Of course, the series expansion (12.5b tacitly assumes that /(±r') possesses con- 
tinuous derivatives of arbitrary order with respect to the Cartesian components of 
its argument ±r' — ±(x',y,z'). Thus, / has to be analytic at ±r', which guaran- 
tees that the series expansion (12.5b converges at least for sufficiently small \r\ > 0. 
Eventual singularities of / determine for which values of the shift vector r the three- 
dimensional Taylor expansion (12.5b converges. 

We could also expand / about r and use ±r' as the shift vector. This would pro- 
duce an addition theorem for f{r±r') in which the roles of r and r' are interchanged. 
Both approaches are mathematically legitimate and equivalent if/ is analytic at r, r', 
and /• ± r' for essentially arbitrary vectors r, / e M^. This is normally not true. Many 
of the functions, that are of interest in the context of molecular electronic structure 
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calculations, have a singularity at the origin. Obvious examples are the Coulomb po- 
tential, which has a pole at the origin, or the commonly used exponentially decaying 
functions as for example Slater-type functions, which have a branch point singularity 
at the origin. 

Accordingly, for the addition theorem of a function /(r ± r'), which is singular for 
r ± r' = 0, the natural variables are not r and r', but the vectors r< and r> satisfying 
r<| = min(r, /) and r>| = max(r, r'). This implies that the expansion formula (12.5b 
should be rewritten as follows: 

f{r^±r:.) ^ £ ^^<'y^^" /(±/->) = e^<-^>/(±r>). (2.6) 

In this way, the convergence of the three-dimensional Taylor expansion is guaranteed 
provided that / is analytic everywhere with the possible exception of the origin. 

From a practical point of view, the translation operator e*"* '^> does not seem to 
be a particularly useful analytical tool. In electronic structure theory, we are usually 
interested in addition theorems of irreducible spherical tensors of the type of (12.11 ). 
which are defined in terms of the spherical polar coordinates r, 9, and (p. Differenti- 
ating such an irreducible spherical tensor with respect to the Cartesian components of 
r> — (x>,3'>,z>) would lead to extremely messy expressions and to difficult techni- 
cal problems. Thus, it is a seemingly obvious conclusion that the translation operator 
gr<v+> Qjjjy provides a formal solution to the problem of separating the variables /•< 
and r> of a function /(/"> ± r>). Nevertheless, this conclusion is wrong. 

The crucial step, which ultimately makes the Taylor expansion method practically 
useful, is the expansion of the translation operator e''<'^> in terms of differential 
operators that are irreducible spherical tensors with a fixed rank ^ G Nq: 

o° ^ ~ 2k ^2k 

e^<-^> ^2nj: J: [^rir<)r ^''"(V>) E 2^^2kkuuk . Z ^'''^ 

i=Om=-e k=0^ i^-(i-/^)(:+k+l 

Here, ^'"(V) is the so-called spherical tensor gradient operator which is obtained 
by replacing the Cartesian components of r = {x,y,z) by the Cartesian components 
of V = {d/dxd/dy,d/dz) in the explicit expression of the solid harmonic ^"(r) = 
r Y^ir/r) which is a homogeneous p olyn omial of degree L More details as well as 



numerous references can be found in 112021 Section 2] 



It seems that the expansion ( I2.7l l was first published by Santos 111731 Eq. (A. 6)], 
who emphasized that this expansion should be useful for the derivation of addition 
theorems, but he appa rently never used it for that purpose. 

In 1 199, 201 ] or in 12021 Section 7] it was shown that two-range addition theorems 
of irreducible spherical tensors are nothing but rearranged three-dimensional Taylor 
expansions about ±r> with shift vector r< . Such an expansion converges pointwise 
and uniformly in the interior of suitable subsets of M^ x M^ whose boundaries are 
defined by the singularities of the function which is to be expanded. If the function 
under consideration has a singularity for r^r' =0, the two-range form of its addi- 
tion theorem is necessary to guarantee pointwise or even uniform convergence in a 
suitable neighborhood of the expansion point. 



11 



There is a different class of addition theorems based on Hilbert space theory. 
Let us assume that {(p'^({i')}nj.m is a complete and orthonormal function set in the 
Hilbert space 

l2(r3) = //: r3^C /'|/(r)|2dV<oc| (2.8) 

of functions that are square integrable with respect to an integration over the whole 
K^ . Since any f EL^ {M? ) can be expanded in terms of the complete and orthonormal 
functions {(pll\;{r)}n,e.m, aone-range addition theoremfor/(r±r') can be formulated 
as follows: 



fir±r') = ^^X/iirOC/W, (2.9a) 

nim 

C':Af;±r') = |[C^,(r)]*/(r±r')dV. (2.9b) 

The expansion ( 12. 9t is a one-range addition theorem since the variables r and r' are 
completely separated in a unique way independent of the lengths of r and r': The 
dependence on r is contained in the functions (p'^^{r), whereas / occurs only in the 
expansion coefficients Cl"f,{f; ±r') which are overlap or convolution-type integrals. 

Since the functions {(p™^(r)}„ £_„, are by assumption complete and orthonormal 
in the Hilbert space L'^(R^), the existence of the one-range addition theorem ( I2.9l l as 
well as its convergence in the mean is guaranteed if / G L^(M^). 

The summation limits in ( |2.9l l depend on the exact definition of the function set 
{Vnii'')}n,(,m- Unless explicitly specified, this article tacitly uses the convention 

oo n-l e 

LLL E ' (2.10) 

nt:m u=l £=0 m= — l: 

which is in agreement with the usual convention for the bound state hydrogen eigen- 
functions. 

It is a typical feature of all expansions of Hilbert space elements in terms of 
function sets, which are complete and orthonormal in this Hilbert space, that they 
do not necessarily converge pointwise but only in the mean with respect to the norm 
of the corresponding Hilbert space. As is well known from the mathematical litera- 
ture, convergence in the mean is weaker than pointwise convergence. Accordingly, 
expansions of that kind are not necessarily suited for a pointwise representation of a 
function. However, as discussed in more detail in Section |6l we must use a weaker 
form of convergence if we want construct one-range addition theorems for functions 
that are singular for r ± r' = 0. Moreover, we do not really need the stronger point- 
wise convergence if we only want to evaluate multicenter integrals with the help of 
addition theorems. 

One-range addition theorems of the kind of (12.91 ) were constructed by Filte r and 
Steinborn |29, Eqs. (5.11) and (5. 12)] and later applied by Kranz andSteinborn 11491] 
and by Trivedi and Steinborn 1 190]. An alternative derivation of these addition theo- 
rems based on Fourier transformation combined with weakly convergent expansions 
of the plane wave exp(±i/? ■ r) with p, r G M^ was presented in il94l Section VII] and 
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in 1113911 . A similar appr oach based on the work of Shibuya and Wulfman 1117611 was 
pursued by Novosadov flSr, Section 3]. 

As discussed in |204, Section 3]), it is also possible to formulate one-range addi- 
tion theorems that converge with respect to the norm of a weighted Hilbert space 

lUm.') = {/:M3^c||w(r)|/(r)|2dV<-}, (2.11) 



where w{r) 7^ 1 is a suitable positive weight function. If we assume that / € L^^,( 
and that the functions {Y'ne{'')}n.l,m ^e complete a nd or thonormal in L^^,(IR^), then 



we obtain the following one-range addition theorem 11204 Eq. (3.6)]): 



/(r±r') = l^ C:/f,w;±r^)w^:,{r), (2.12a) 



nim 



C"„>(/,w;±r') = J [v/;;,(r)]*w(r)/(r±r')dV. (2.12b) 

A one-range addition theorem of the type of either ( |2.9l l or (12.121 1 for a function 
/ : M^ ^ C is a mapping R^ x M^ — > C. Compared to the better known two-range 
addition theorems like the Laplace expansion ( 12.2b . which depend on r and r' only 
indirectly via r< and r>, one-range addition theorems have the highly advantageous 
feature that they provide unique infinite series representations of functions /(r ± r') 
with separated variables r and r' that are valid for the whole argument set M? xM?. 

In his fc-dependent one-range addition theorems for Slater-type functions and re- 
lated functions, Guseinov used as a complete and orthonormal set the following class 
of functions 1461 Eq. (1)], which - if the mathematical notation (ID. 2b for the general- 
ized Laguerre polynomials is used - can be expressed as follows 1.204. Eq. (4.16)]: 



kXiir.r) = 



(27)^+3(„_^_l)!]'/2 



{n+e+k+iy. 



.-r'-r(2m+2)^27r)^™(27r), 



^n-l-\ 



neN, n>^+l, A:=-l,0,l,2,..., 7>0. (2.13) 



The restriction to integral values of A; = — 1,0,1,2,... is unnecessary. The mathemat- 
ical definition ( ID. 2b of the generalized Laguerre polynomials f/„ (z) permits nonin- 
tegral superscripts a. Thus, the condition fc = — 1 , 0, 1 , 2, . . . in (12.13b can be replaced 
by A: £ [— 1 , 00) . If this is done, one only has to replace the factorial (n + ^ + A: + 1 ) ! in 
(12.13b by the gamma function r{n+ l + k + 2). 

These functions are orthonormal with respect to the weight function w(r) = r^ 
(compare also ||46l Eq. (4)]): 

/ [^"0(7, r)] * r'^ ,f„?; (7, /■) dV = 5„„, 5u' 5w . (2. 14) 

Accordingly, these functions are complete and orthonormal in the weighted Hilbert 
space 

i'.(IR') = {/:K'^c||/|/(r)|2dV<-}, ^=-1,0,1,2,.... (2.15) 
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As discussed in the text following 11204 Eq. (4.20)], the functions k^ll\{'Y^r) can 



- depending on the value of the free parameter A: = — 1,0,1,2,... - reproduce several 
other physically relevant complete and orthonormal function sets. 

If we set A; = in ( 12.151 ). we retrieve the Hilbert space L^ {M? ) of square integrable 
functions defined by ( 12.81 1 which because of the Born interpretation of bound state 
wave functions is the most natural choice for the representation of effective one- 
particle wave functions in electronic structure calculations. 

For k ^ 0, there is, however, a problem since we neither have L?{M?) C L\{M?) 
nor L^^t{R^) C L^{R^). Thus, the Hilbert spaces L^{R^) and L^^(M^) are for /fe 7^ 
inequivalent. This applies also to approximation processes which (should) converge 
with respect to the norms of these Hilbert spaces (see for example [206, Appendix 
E]). 

Conceptually, the derivation of a one-range addition is a triviality. This follows 
from the obvious fact that an arbitrary Hilbert space element can be expanded in 
terms of a function set that is complete and orthonormal in this Hilbert space. 

Unfortunately, this does not imply that the derivation of a one-range addition 
theorem is necessarily a simple task. The challenging part is the construction of 
computationally convenient expressions for the overlap integrals ( I2.9bb or (I2.12bl ). 
respectively. In realistic applications, we cannot not tacitly assume that the use of 
one-range addition theorems in multicenter integrals necessarily leads to rapidly con- 
vergent expansions (compare for instance the convergence rates reported by Trivedi 
and Steinborn |190]). Therefore, we must be able to compute the overlap integrals 
( I2.9bb or ( I2.12bl ) both efficiently and reliably even for possibly very large indices. 

Let me conclude this Section with some short remarks on the relative advantages 
and disadvantages of some of the various sets of exponentially decaying function sets 
that are described in the literature. 

The best known and most often used exponentia lly decayi ng functions are un 



doubtedly Slater-type functions ( 11.11 ) introduced in 111781 117911 . Slater himself had 
argued that sufficiently accurate approximations to numerically determined solutions 
of effective one-particle Schrodinger equations can be obtained even if the radial 
nodes of these solutions are completely ignored 11781 p. 57]. Slater's prime concern 
was not accuracy but analytic simplicity fl79, p. 42]. In view of the limited compu- 
tational resources at that time. Slater's pragmatic attitude certainly made sense. 

Because of their remarkably simple structure. Slater-type functions are often con- 
sidered to be the most basic prototypes of all exponentially decaying functions. How- 
ever, the simplicity of Slater-type functions in the coordinate representation is decep- 
tive. Slater-type functions have been used with considerable success in the case of 
atomic electronic structure calculations, where multicenter integrals do not occur. 

But there is considerable evidence that in the case of multicenter integrals it is 
the simplicity of a functio n in the momentum representation that really matters. As 
for example discussed in f216T, the Fourier transform of a Slater-type function is a 
comparatively complicated object having the same level of complexity as the Fourier 
transform of a bound state hydrogen eigenfunction (see for example 11941 Section 
IV]). Therefore, it is certainly worth while to look for alternative exponentially de- 
caying function sets with simpler Fourier transforms and also more convenient prop- 
erties in multicenter problems. 
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Inspired by previous work of Shavitt il75l Eq. (55) on p. 15], Steinborn and Filter 



1 1821 Eqs. (3.1) and (3.2)] introduced the so-called reduced Bessel function 

icyiz) = (l/Tty^h^Kyiz). (2.16) 



Here, Ky{z) is a modified Bessel function of the second kind 11154 p. 66]. If the 
order v is half-integral, v = m + 1/2 with m G No, the reduced Bessel function is an 
exponential multiplied by a terminating confluent hypergeometric series i^i (see for 
example |l217l Eq. (3.7)]): 

L+i/iiz) =2"'{l/2Ue-'iFi{-m;-2m;2z). (2.17) 



As discussed in more detail in Section |6] or in 12081 Section 4], Steinborn and Filter 
became interested in reduced Bessel functions because of a known Gegenbauer-type 
addition theorem which allowed a simple derivation a two-range addition theorem for 
reduced Bessel functions w ith h alf-integral orders (see for example |182, Eq. (3.4)] 
or, as an improved version 12181 Eq. (5.5)]). 

In connection with convolution and Coulomb integrals. Filter and Steinborn later 
introduced the so called B functions as an anisotropic generalization of the reduced 
Bessel functions with half-integral orders 1281 Eq. (2.14)]: 

fi™,(/3,r) = [2"+^(« + ^)!]-i^„_i/2(j3r)^"'(j3r), neZ. (2.18) 

B functions are fairly complicated mathematical objects, and ( 12.181 ) and ( 12.171 ) 
imply that a B function can be expressed as a linear combination of Slater-type func- 
tions with integral principal quantum numbers. Hence, it is not at all clear why the 
use of the comparatively complicated B functions should offer any advantages over 
Slater-type functions which possess an exceptionally simple explicit expression in 
the coordinate representation. 

Let us for the moment assume that we form some finite linear combinations of 
Slater-type functions and that we do some mathematical manipulations with this lin- 
ear combination. Normally, the complexity of the resulting expression increases, de- 
pending on the number of Slater-type functions occurring in the linear combination. 
In fortunate cases, however, it may happen that most terms of the resulting expression 
cancel exactly. Thus, a significant reduction of complexity is also possible if we form 
appropriate linear combinations. 

The Fourier transform of a B function seems to be such a fortunate case since it 
is of exceptional simplicity among exponentially decaying functions: 

BUa,p) = (27r)-3/2 / e-'P-'-B™,(a,r)dV 



a 

\a 



2n+e-l 



(2/^)'^ \.2Tn2i»+.+i ^/"(-^P)- ^2.19) 



This is the most consequential and also the most often cited result of my PhD thesis 



111931 Eq. (7.1-6) on p. 160]. Later, dTTgl l was pubUshedin [216, Eq. (3.7)]. Indepen- 
dently and almost simultaneously, ( 12.19b was also derived by Niukkanen 1.156k Eqs. 
(57) - (58)]. 
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The exceptionally simple Fourier transform (12.19b explains why multicenter inte- 
grals of B functions are often simpler than the corresponding integrals of other expo- 
nentially decaying functions (see for example 1 170, 171, 208] and refe rences ther ein). 



It also explains why it was comparatively easy to derive two-range 112011 121811 and 
one-range 12911 addition theorems forB functions. 

The exceptionally simple Fourier transform ( 12.191 ) also explains why other expo- 
nentially decaying functions can be expressed in terms of B functions. For example, 
a Slater-type function with an integral principal quantum number can be expressed 
by the following finite sum of B functions 1281 Eqs. (3.3) and (3.4)]: 

Z„">(/3,r) = 2" E ^_,) A~[n^i-m.i~[n^mU 

(7>0 '^■ 

x(«-(T)!B;r_,_„/j3,r). (2.20) 

If n — ^ is even, the Pochhammer symbol ( — [n — i]/2)a causes a termination of the 
(7 summation after a finite number of steps, and if n — ^ is odd, this is accomplished 
by the Pochhammer symbol (— [n — £ — l]/2)<j. 

It is also possible to express Guseinov's complete and orthonormal functions 
j(,f^™,(7, r) defined by ( 12.131 1. which play a central role in Guseinov's work on one 



range addition theorems, as a finite sum of B functions. We only have to use 11193 . 
Eq. (3.3-35) on p. 45] 

to obtain 



f+^{n+e+k+iy.Y^^ {2n+k+i)r{i/2){i+iy. 



2*^+1 („-£_i)! / r{e+2+k/2)r{e+[k+5]/2) 



n-e-\ 



h vl{l + 2 + k/2\{l+[k + 5]/2)^ 5v+i.(7,r). (2.22) 

Consequently, explicit expressions for multicenter integrals and addition theo- 
rems of most exponentially decaying functions - among them Slater-type functions 
with integral principal numbers and Guseinov's functions - can be derived via the 
often simpler analogous results for B functions. 



3 Guseinov's Rearrangements of One-Range Addition Theorems 



In ll4ll - l44ll46ll . Guseinov derived one -range addition theorems for Slater-type func- 
tions XNLiP J *" i '■') with integral or nonintegral principal quantum numbers A^ by ex- 
panding them in terms of his complete and orthonormal functions {i^f^™ (7, '')} « 
defined by ( 12.131 1 with an in general different scaling parameter 7 ^ j3 > 0: 

x}f4P,r±r') = EXf;f(7.^.±'-')A.'0(7,/-), (3.1a) 

nim 

Xffir^P:^'-') = [[,%y{y,r)]*r'x^Afi.r±r')Ah. (3.1b) 
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As long as the principal quantum number A^ is not too negative, which will be tacitly 
assumed in the following text, the Slater-type function ;i;^^(j3,r ± r') belongs to the 
weighted Hilbert spaces L\(M?) defined by ( 12.15b with fc = -1,0, 1,2,.... Accord- 
ingly, the ^-dependent one-range addition theorems ( 13.11) . which are special cases of 
the general addition theorem ( 12.121 ). converge in the mean with respect to the norms 
of their corresponding weighted Hilbert spaces L\ (R^)- 

As akeady remarked at the end of Section |2] the central computational problem 
occurring in the context of Guseinov's one-range addition theorems ( 13.11 ) is the effi- 
cient and reliable evaluation of the overlap integrals (I3.1bb . 

Normally, it is much easier to compute overlap integrals of exponentially decay- 
ing functions with equal scaling parameters j3 = 7 than with different parameters 
j3 ^ 7. However, Guseinov wanted to have this additional degree of freedom for cer- 
tain applications. For example, in |66|] Guseinov represented the Coulomb potential 
1 / 1 r — r' I as the limiting case )3 ^- of the Yukawa potential exp (— )3|r — r'|)/|r — r'|, 
which is proportional to the Slater-type function Xoo{P,r — r').ln this way, Guseinov 
ll66ll formally obtained a one-range addition theorem of the Coulomb potential in 
terms of his functions i.f„™ (7, /"). But this approach leads to other convergence and 
existence problems [204, 206] which will be discussed in more details in [211]. 

For the evaluation of the overlap integrals ( I3.1bb occurring in his addition the- 
orems ( 13. lb ). Guseinov uses a simple approach which was first used by Smeyers 
LI 801 in 1966. To the best of my knowledge, this approach was adopted in 1978 by 
Guseinov 114 li Eqs. (6) - (8)] and consistently used in his later publications. Unfortu- 
nately, Smeyers' simple approach is not necessarily good, as akeady emphasized in 
1982 by Trivedi and Steinborn [190, pp. 1 16 - 1 17]. 

It follows at once from the explicit expression ( ID.2b of the generalized Laguerre 
polynomials that Guseinov's functions can be expressed as finite sums of Slater-type 
functions with integral principal quantum numbers: 



k%":ei7,r)= t kGtTiiT)x':Ar,r), (3.2a) 

v=e+i 



,Gf\y) = 2^ 



{2yf+\n + (. + k+\)\ 



(«-£-!)! 






(3.2b) 



Accordingly, the overlap integrals (I3.1bb can be expressed as finite sums of overlap 
integrals of Slater-type functions: 

X:t^ir,li,±^) = r' t A"lU7)S:t,l„i7,li,±r'), (3.3a) 

S^J'^{r,l3,±r') = f[x':/Y,r)Yx^Ali,r±r')d\. (3.3b) 



Inserting this into the addition theorems (13.1b yields (compare 11461 Eqs. (15) and 
(16)]): 

z^j^,/-±/) = 7-'^E.f^;"(7,'-) t A"-tiir)s^vtuay^P^±'''). (3.4) 

nim v=i+l 
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Guseinov's approach seems to have the advantage that existing programs for over- 
lap integrals of Slater-type functions can be used for the evaluation of the overlap 
integrals ( I3.1bb involving Guseinov's functions. Unfortunately, this seemingly con- 
venient approach can easily lead to stability problems. Numerical instabilities are 
extremely likely in expressions based on (13.21 ) if the indices n of the Guseinov func- 
tions 1;.^^™, (7, r) are large. 

In his desire to reduce his whole formalism of one-range addition theorems to 
Slater-type functions with integral principal quantum numbers, Guseinov even ex- 
pressed the functions ^^f^'" (7, /■) on the right-hand side of ( 13.41) by Slater-type func- 
tions according to ( 13.2b (compare i4a Eqs. (14) - (16)]): 

X^J^,r±r') = 7-'L I M'-li(7)Zv5i(7,'-) 

nim v'=l+l 

X t kG^:-Uy)S':'^,%{r^P,±r'). (3.5) 

In this version of Guseinov's A:-dependent addition theorems, the complete and or- 
thogonal functions ^,f^J",(7, r) do not occur any more and are replaced by Slater-type 
functions X'^^iY-'') with integral principal quantum numbers n > £+1. Apparently, 
Guseinov considered this to be a major achievement. 

But there is a price to be paid. Firstly, introducing one finite inner sum after the 
other does not look like a promising computational strategy. Secondly, the numerical 
stability of Guseinov's transformations is questionable, in particular in the case of 
(very) large quantum numbers. 

Whenever we do calculations of nontrivial complexity with a fixed precision, 
numerical stability is invariably a very important issue. Let us for example assume 
that we work with a series expansion whose coefficients are given by complicated 
expressions such as multiple nested finite sums (this characterization certainly applies 
to Guseinov's rearranged addition theorems). If alternating signs occur in these nested 
sums, a loss of significant digits cannot be ruled out, and in the case of (very) large 
indices, a catastrophic accumulation of rounding errors may even produce completely 
nonsensical results. 

Guseinov's rearranged addition theorems are based on the simple fact that Gu- 
seinov's orthogonal functions ^(,^^'",(7, r) can according to (13.2b be expressed as finite 
sums of Slater-type functions with integral principal quantum numbers. The linear 
combination ( 13.21 ) is nothing but the explicit expression ( ID.2b of a generalized La- 
guerre polynomial in disguise. The coefficients in (ID. 21) and thus also the coefficients 
in ( 13.2b have strictly alternating signs. 

Sign alternation of the coefficients of the classical orthogonal polynomials is a 
necessary requirement for an orthogonality relationship of the type of ( ID. 3b . But this 
sign alternation can easily lead to a catastrophic accumulation of rounding errors, if 
we try to compute an orthogonal polynomial with a large index from its explicit ex- 
pression. In practice, one usually avoids explicit expressions. It is much better to eval 



uate orthogonal polynomials via their three-term recursions (see for example 1163 , 



Eq. (18.9.1) and Table 18.9.1]). For example, Gautschi, who is generally considered 
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to be one of the leading experts on the stabiUty of Hnear recurrences, wrote in |33, p. 
277]: 

It is our experience, and the experience of many others, that the basic three- 
term recurrence relation for orthogonal polynomials is generally an excellent 
means of computing these polynomials, both within the interval of orthogo- 
nality and outside of it. 

A FORTRAN SUBROUTINE OTHPL, which evaluates classical orthogonal poly- 
nomials via their three-term recursions 1 163, Eq. (18.9.1) and Table 18.9.1], is listed 
in 1I220I pp. 23 - 24]. Another important aspect is that a recursive evaluation is much 
more efficient, in particular if whole strings of orthogonal polynomials have to be 
evaluated simultaneously. 

As long as the opposite is not expUcitly proven, it makes sense to be cautious and 
to assume that an expression like (I3.2l i. whose coefficients ^.Gy (7) have strictly al- 
ternating signs, inherits the stability problems of the explicit expression ( ID.2l i of the 
generalized Laguerre polynomials from which it was derived. Similarly, the coeffi- 
cients iS^%'_f;_\{y) and ^G^_i,_y{y) in the inner v' and v sums in ( I3.5l l have strictly 
alternating signs. 

1 am not at all convinced that inner sums based on ( ID.2l i. which occur as coeffi- 
cients in infinite series expansions like the one on the right-hand side of ( 13.51 ) or later 
in (13.71 ) and ( 13.8b . can be computed in a numerically stable way for large values of the 
outer summation indices n and I. This has to be investigated (much) more thoroughly. 

Another problem is that conventional programs for overlap integrals of Slater- 
type functions, as they are for instance used in semiempirical calculations, normally 
cannot be used in the case of (very) large principal and angular momentum quantum 
numbers. Therefore, one should look for alternative expressions for the overlap inte- 
grals ( I3.1bl ). which are not based on (I3.2l l and which permit an efficient and reliable 
evaluation even for (very) large values of the indices n and (,. 

Slater-type functions are complete but not orthogonal. But this nonorthogonal- 
ity can also cause problems. There is a practically very consequential aspect of or- 
thogonal expansions which is often not sufficiently appreciated. Let us assume that 
/ belongs to some Hilbert space ,^ with inner product (1) and that the functions 
{<Pn}Ti=Q ^^ complete and orthonormal in ^. Then, as discussed in Appendix ICl / 
possesses the orthogonal expansion ( IC.4b which converges in the mean with respect 
to the norm || • || of Jf . Moreover, the expansion coefficients {<Pn\f) satisfy Parseval's 
equality (see for example 1 166, Eq. (11.2) on p. 45]) 

\\f\? =JL\{<Pn\f)\^- (3.6) 

Parseval's equality implies that the inner products {<Pn\f) are bounded in magnitude 
and that they vanish as « —5- 00. This may well be the main reason why orthogonal 
expansions tend to be computationally well behaved. 

In the case of nonorthogonal expansions of the type of ( IC.3b , quite a few com- 
plications can happen. We cannot tacitly assume that the expansion coefficients C„ in 
(IC.3b are necessarily bounded in magnitude. These coefficients can have alternating 
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signs and even increase in magnitude with increasing index n (examples can for in- 
stance be found in 1144 Table I on p. 166] or fir, Appendix E on pp. 162 - 164]). 
Such a behavior of the expansion coefficients can easily lead to a cancellation of 
significant digits or even to a catastrophic accumulation of rounding errors. 

Guseinov approximated the Slater-type function XnlH^ , ''i J^) by a truncation of 
the addition theorem ( 13.5b including only the first ^ terms of the outer n summation 
and defined the complete addition theorem as the limit yy -^ oo of his ^/T-dependent 
truncation (compare |46, Eq. (13)]): 

Z^,(i3,r±r') = 7-' lim ft I 

X t M'li(7)Zv"5,(7,'-) t .Gi'^li(7)S:;^:^,„.(7,^,±'-'). (3.7) 
v'=e+i v=i+i 

In the .yK-dependent part of this expression, Guseinov changed the order of summa- 
tions and formally expressed this finite sum as a linear combination terms of Slater- 
type functions with integral principal quantum numbers |46, Eq. (15)]. In this context, 
it is in my opinion advantageous to change the order of the finite n and £ summations 
according to ^;;ti L"ld ^ L'Cq^ l^n=e.+\- Then, we obtain: 

Z^J^,r±r') = y->^ Jim 'f ^ I ZS+i(7,/-) 



■A 



X 



l=a m=-('t=li+\ 



E.G^':fi(7) E .G^.^(7)S^:^j7,/3,±r'). (3.8) 



As long as ,yV is finite, the transformation of the c/K-dependent part of (I3.7l i to the 
^-dependent part of (I3.8l l is legitimate. But it is not guaranteed that this transforma- 
tion is still legitimate in the limit .jV — s> oo and that one-range addition theorems can 
be reformulated as expansions in terms of nonorthogonal Slater-type functions with 
integral principal quantum numbers. Further details will be given in Section|6] 

One-range addition theorems for Slater-type functions are fairly complicated math- 
ematical objects, whose series coefficients are essentially overlap integrals. Thus, a 
detailed analysis of the existence and convergence properties of such an addition 
theorem is certainly a very demanding task. Consequently, it makes sense to look 
for simplifications and to pursue an indirect approach. The situation becomes much 
more transparent if we do not look at addition theorems for Slater-type functions 
XiJii^ , r ± r'), but rather at their angular projections 

S,^f'"(^,r,±r') = / [y/"(rA)]*;t^^(i3,r±r')dV. (3.9) 

J\r\ = \ 

If we insert the addition theorem ( 13.11 ) into this integral and perform the integration 
over the surface of the unit sphere in M?, we obtain with the help of the orthonormal- 
ity of the spherical harmonics an expansion in terms of the radial parts of Guseinov's 
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functions ,5.^^™ (7? '")> which can be reformulated as an expansion in terms of general- 
ized Laguerre polynomials: 



yr^N,L.Mio .. i _/^ X"' vN.LM 

e Li/j 



(i3,r,±r') = l^ ,x::;^;f (7,/3,±r') 

n=l+\ 



{2yf+\n~e.^l)\'^^'^ 



{n + i + k+\)\ 



(yrfL^^lt-fhlYr). (3.10) 



A further substantial simplifications takes place if we do not consider the fairly com- 
plicated one-range addition theorems but rather their much simpler one-center limits 
r' = 0. Then, for fixed j3 , 7 > the overlap integrals ^X^^ '^ ^'^ (7) j3 , ^r') defined by 
(I3.1bb simplify to become numbers. Moreover, because of the orthonormality of the 
spherical harmonics only the overlap integrals with (i^L and m = M can be nonzero. 



4 The Transformation of Laguerre Series to Power Series 

If Guseinov's rearrangements of the one-range addition theorems (13.1b for Slater- 
type functions Xnl iP^''^''') ^"^ legitimate, we ultimately obtain an expansion of a 
Slater-type function with an in general nonintegral principal quantum number N eM.\ 
N in terms of Slater-type functions {XneiY^'')} / with integral principal quantum 
numbers n and a common scaling parameter 7 > 0: 

x}f,LiP,r±r') = ^ xS'^)(i3,7,±r')z„":,(7,r). (4.1) 

If we set r' = 0, this expansion simplifies and we obtain ( 11.2b . 

If the expansion ( 14.11 ) exists and converges for the whole argument set M? x R^, 
then it is obviously a one-range addition theorem. Moreover, ( 14.11 ) can be reformu- 
lated as a power series in yt" by shifting the exponential part exp(— 7r) of the Slater- 
type functions Xnei'Y''') '^° ^^^ left-hand side: 

ey'x^AP,r±r') ^ ^ xi^,':;^\[i,r,±r'){yrr-'Yr{r/r). (4.2) 

This expansion can be transformed to an infinite number of ^-dependent power series 
in yr for the angular projections S^ ^^ ' (j3 , r, ±r') defined by ( 13.91 ): 

er^<;f'^(i3,r,±r') = f xiy/'(j3,7,±r') (7'-r ' ■ (4.3) 

A comparison of ( 13.101 ) and ( 14.3b shows that Guseinov's rearrangements ulti- 
mately corresponds to the transformation of a Laguerre series 

f{z) = f:Ai")4"'(z), (4.4a) 
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which converges in the mean with respect to the norm of the weighted Hilbert space 
L^-z^a ([0,°°)) defined by ( ID.8I 1, to a power series 

fiz) - E7«^", (4.5) 

n=0 

which - if it exists - converges pointwise in a suitable subset of the complex plane. 

There is, however, the problem that analyticity in the sense of complex analy- 
sis and the existence of a Laguerre expansion are completely unrelated concepts. 
Therefore, we cannot expect that the transformation from (14.41 ) to ( 14.51 ) is necessarily 
possible in the case of an essentially arbitrary function / : C — > C. 

The convergence and e xiste nce problems, which occur in this context, were re- 



cently studied in depth in 1120711 . There, some simple sufficient conditions based on 
the decay rate and the sign pattern of the Laguerre series coefficients A„ ^ as « ^ 00 
were formulated. These conditions, which extend previous results by Gottlieb and 
Orszag |35, p. 42] and by Doha f25l, p. 5452], respectively, make it possible to decide 
whether the transformation of a Laguerre series of the type of ( 14.41 ) to a power series 
produces a mathematically meaningful result or not. 

In order to understand better the possible pitfalls and dangers of such a transfor- 



mation, let us first consider a partial sum of the infinite Laguerre series (14.4b : 

M 

/m(z) = L x!"'>Li"\z) (4.6) 

In this partial sum, the power z'' with < p <M occurs in the Laguerre polynomials 
Lp (z), ip4_i (z), . . . , Ll^ (z). Thus, we have add up all contributions with v = p on 
the ri ght-hand side of (14.7b to obtain the coefficients of z''. A short calculation yields 
ll20lEq,(3.13)]: 

^ ^ i-zY y (g + v + l), („) 

v=0 ^- li=0 f^- 

Since the truncated Laguerre series /m(z) is simply a polynomial of degree M in 
z, it is always possible to reformulate it by interchanging the order of the nested finite 
sums. No convergence and/or existence problems can occur. 

Unfortunately, this does not mean that the transformation from (14.6b to (14.8b is 
also possible in the limit M ^> 00. In this case, the finite inner sum on the right- 
hand side of ( 14.8b becomes an infinite series which can a diverge if the coefficients 

(oc) 
X„ do not decay sufficiently rapidly as n — > 00. Moreover, the fact that /m(z) is a 

mathematically meaningful object does not imply that /(z) possesses a power series 

expansion of the type of ( 14.5b which converges in a suitable subset of the complex 

plane. The expression (14.8b for /m (z) possesses the following general structure: 

M 

/m(z)= l^rl^hr (4.9) 

M=0 
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The coefficients 7,, depend explicit on the summation limit M, i.e., we have in 
general 7„ ^ 7,, y^Yn i^ ■■■ for all n e Nq. Thus, we have to show explicitly 
that the 7, ' possess for all « £ No well defined limits 7„ = y^ — limM^oo 7, . 
Only if all these limits exist, /(z) can possess a convergent power series expansion 

about z = 0. This, however, requires that the Laguerre series coefficients A,, decay 
sufficiently rapidly as « — )► 00. 

Let us now consider the transformation of the infinite Laguerre expansion (14.4b . 
The power z'^ with p >0 occurs in all Laguerre polynomials L„ (z) with n> p. We 
obtain the corresponding power series coefficient Yp, if we add up all contributions 
with V — pin the infinite series 



A short calculation yields |I207I Eq, (3.14)]: 

/w-i;^E^^^i^<.i- (4.11) 

v=0 ^- ^=0 f^- 

If we compare ( I4.8l l and ( I4.11l i. which corresponds to the limit M -^ 00 in ( 14.81 ), 
we immediately see that the rearrangement of an infinite Laguerre expansion is not 
necessarily possible since the power series coefficients 7, are now given by infinite 
series. If the inner ji series in ( 14.11b do not produce finite results, we end up with a 
formal power series with expansion coefficients that are infinite in magnitude. 

Thus, the existence of the power series (14.5b for /(z) depends crucially on the 
decay rate and the sign pattern of the Laguerre series coefficients X,\ as n — s> 00. In 



1 20711 . three different prototypes of large index behavior of the coefficients X„ were 



discussed: 

- algebraic decay with ultimately monotone signs, 

- exponential or factorial decay, 

- algebraic decay with ultimately strictly alternating signs. 

If the coefficients X„ decay algebraically as n — > 00 and ultimately have mono- 
tone signs, the inner ji series in (14.111 ) diverge for sufficiently large outer indices v 



1 2071 Section 4]. Thus, a function /(z) represented by a Laguerre expansion with 



1 eel 
ultimately monotone and algebraically decaying series coefficients A„ cannot be 

analytic in a neighborhood of the origin z = 0. 

If the coefficients A„ decay exponentially or even factorially, the inner jj. series 
in (14.11b converge (often quite rapidly), and the function/ (z) under consideration is 
analytic in a suitable neighborhood of the origin z = (ZCQ Section 5]. This is pretty 
much the most advantageous situation we can hope for 

These conclusions are in agreement with short remarks by Gottlieb and Orszag 
ySi p. 42] and by Doha [25, p. 5452], respectively, who had stated without detailed 
proof that such a Laguerre series converges faster than algebraically if the function 
under consideration is analytic at the origin. However, both Gottlieb and Orszag 11351 
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p. 42] and Doha 11251 p. 5452], respectively, had failed to see that there is a different 
scenario, which also produces convergent power series expansions about z = 0. 

(ex) 
Let us assume that the coefficients A„ decay algebraically in magnitude and ul- 
timately have strictly alternating signs. Then, the inner jj. series in ( 14.111 ) still diverge 
at least for sufficiently large values of the outer index v. However, this divergence 
does not imply that a convergent power series about z = does not exist. The strictly 
alternating signs make it possible to associate something finite to the divergent ji 
series by employing suitable summation techniques for divergent series. A highly 
condensed review of summation techniques with a special emphasis on sequence 
transformations can be found in |207, Appendices A and B], and in |209] there is 
an admittedly incomplete discussion of the usefulness of sequence transformations 
beyond the summation of divergent series. 



In 112071 Section 4], several known generating functions of the generalized La- 
guerre polynomials were recovered by summing divergent inner jj. series representing 
power series coefficients. In these comparatively simple cases, the divergent inner /i 
series could always be expressed by generalized hypergeometric series with argument 
— 1 . However, these generalized hypergeometric series only converge in the interior 
of the unit circle. Thus, analytic continuation from the interior of the circle of conver- 
gence of such a hypergeometric series to its boundary was sufficient to accomplish a 
summation. In the case of simpler hypergeometric series, this is not too difficult. 

A skeptical reader might therefore conclude that a summation is on ly feasible in 
the case of comparatively simple problems. This is not true. In flOr, Section 6] it 
was shown that it is also possible to sum divergent inner fi series by purely numerical 
techniques. Particularly good results were obtained by the nonlinear S^ transforma- 
tion which I had introduced in Iil95t Eg. (8 .4-4)] (a highly condensed review of the 
historical development was given in ll210i Section 2]). Some authors c all this J>^ 
transformation the Weniger transformation (see for example lll3Ml5[ l22l Il52l 118811 



or 1134 Eq. (9.53) on p. 287]). This terminology w as als o used in the recently pub- 
Ushed NIST Handbook of Mathematical Functions il63l Chapter 3.9(v) Levin's and 
Weniger's Transformations] (see also the companion NIST Digital Library of Math- 
ematical Functions under http : //dlmf . nist . gov/3 . 9#v). 

The example of divergent but summable inner jj. series shows that fairly sophisti- 
cated mathematical techniques may be needed to accomplish the transformation of a 
Laguerre expansion to a power series. 

It should be noted that the structure of the formulas, which effect the transforma- 
tion of a Laguerre series to power series, is of a more general nature, and that very 
similar transformation formulas occur also in completely different contexts. In 121(1 
Appendix B] it was shown that formulas having analogous structures occur as long 
as the transformation matrices are triangular and satisfy certain orthogonality con- 
ditions. In 1:210], transformations between factorial series and inverse power series 
were considered, and the transformation matrices involved Stirling numbers of the 
first and second kind, respectively. 
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5 One-Center Expansions for Slater-type Functions 

As already mentioned above, the one-center expansion for Slater-type functions with 
in ge neral nonintegral principal quantum numbers used by Guseinov and Mamedov 



1 1221 Eq. (4)] was originally derived by Guseinov |46, Eq. (21)] as the one-center 
Umit of a rearranged addition theorem for Slater-type functions of the type of ( I3.8l l 
146, Eq. (15)]. It is, however, both simpler and for our purposes much more in- 
structive to derive first one-center expansions of Slater-type function with both in- 
tegral and nonintegral principal quantum numbers in terms of Guseinov's complete 
and orthonormal Laguerre-type functions {k^^'dY,!")} » defined by ( 12.131 ). In the 
second step, it is investigated whether and under which conditions these orthogo- 
nal expansions can be transformed to expansions in terms of Slater-type functions 
{XneiYi'')} I with integral principal quantum numbers n. 

A convenient starting point for the construction of the one-center expansion of a 
Slater-type function with an essentially arbitrary princip al qua ntum number in terms 
of Guseinov's functions is the following Laguerre series 12041 Eq. (6.12)]: 



■p»«z = n „,,^-«-p 



zPe"^ = (1-m) 



_i r(a + p + l) 
r(a + i) 



X £2^1 (-«,a + p + l;a + l;-^j4"'(z), 

peM\No, 5R(a + 2p)>-l, mg (-00,1/2). (5.1) 

The restriction u E {—°°, 1/2) is necessary to guarantee the existence of some inte- 
grals occurring in the derivation of this expansion. If we assume p = m with m £ No, 
no immediately obvious simplification occurs in ( 15.1b . But for m = 0, we obtain a 
terminating Gaussian hypergeometric series 2^1 with unit argument that can be ex- 
pressed in closed form with the help of Gauss' summation theorem [154, p. 40] and 
we obtain a much simpler expansion |26, Eq. (16) on p. 214]; 



zP = 



^^^Eilj^^'"'fe)- P-«»> ^(a+2p,>->.(5.2, 



If we now set p = m with m G No, the infinite series terminates because of the 
Pochhammer symbol (— m)„ and we obtain instead a finite sum: 

z"' = (a + l),„^ > ;" L,r'(z), meNo, 3iia) + 2m>-l. (5.3) 

„=0 [(^+ '-Jn 

The Laguerre expansions listed above can all be transformed to expansions of 
Slater-type functions in terms of Guseinov's functions. The expansion ( 15. Il l yields 
the most general case, i.e., we obtain a A^-dependent one-center expansion of a Slater- 
type function XNiiP ' *") ^i'^h a nonintegral principal quantum number A^ G M \ N in 
terms of Guseinov's functions with different scaUng parameters )3 7^ 7 > 12061 Eq. 
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(5.9)]: 



[j3 + Y]N+L+k+2 {2L + k + 2y. 



v=0 



{v+2L+k+2y: 



1/2 



fv'lL+lx(7,A-) 



X2F1 (-V,N + L + k + 2;2L + k + 3-J^] . (5.4) 

V P + rJ 

This expansion can also be derived by performing the one-center limit r' = in the 
addition theorem ( 13.1b . 

If we assume in (15.4b that the principal quantum number A^ of the Slater-type func- 
tion is a positive integer satisfying A^ > L+ 1, no immediately obvious simplification 
occurs. But in the case of equal scaling parameters j3 = 7 > 0, whic h cor responds to 
M = in ( 15.1b . the expansion ( 15.4b simplifies considerably, yielding 112061 Eq. (5.7)]: 

XNA^.r) = ^ ^^^^_^ r[N + L + k + 2) 



X 



(-A^ + L+l)v 



'^v+L+1,l(^''")' 



i-i r/ N nl/2 k^V+L+\.L\ 

v=o [(v + 2L + /fc + 2)!v!] ' 



A?eM\N, j3>0, /t=-l,0,l,2,.... (5.5) 

If A^ G N and A^ > L + 1 , the infinite series on the right-hand side terminates because of 
the Pochhammer symbol {—N + L+ l)v. The resulting finite sum can also be derived 
directly via ( 15.31 ). 

With the help of the sufficient conditions formulated in 120711 and discussed in 
SectionlH we can analyze whether the Laguerre-type functions j^'f^'"^'(7, r) in the one- 
center expansions given above can be replaced by Slater-type functions X'nciYi'') 
with integral principal quantum numbers. But it is simpler and also more transparent 
to investigate instead whether the equivalent Laguerre expansions (15.1b . ( 15.21 ). and 
( 15.3b can be transformed to power series. 

The most simple situation occurs if we replace in the finite sum ( 15.31 ) for the 
integral power z"' with m € No the generalized La guerr e polynomials by powers via 



( ID.2I ) and rearrange the order of the summations 12071 Eqs. (3.5) and (3.6)]. Then 
we arrive at the trivial identity z"' = z"\ which proves the correctness of ( 15.31 ) but 
provides no new insight. 

In the case of the Laguerre expansion ( 15.2b for a nonintegral power z^ with p G 
K \ No, we face a completely different situation. Firstly, we have no a priori reason to 
assume that (15.2b converges pointwise. Hilbert space theory only guarantees that this 
expansion converges in the mean with respect to the norm of the weighted Hilbert 
space L^_j „([0,oo)). Moreover, the power function z^ is not analytic at z = in 
the case of nonintegral p G M \ No. Thus, a mathematically meaningful power series 
about z = cannot exist, and the transformation of the Laguerre expansion ( 15.2b via 
(14.11b also must fail. 
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This is indeed the case. We only need the asymptotic approximation yj, Eq. 
(6.1.47) on p. 257] 

r{z + a)/r{z + b) ^ z"-''[l+0{l/z)], z-^-o, (5.6) 

to obtain the following asymptotic estimate for the coefficients in ( I5.21 l: 

r(a + i) (a + i)„ r(-p) L V ;j ' 

This asymptotic estimate shows that the coefficients in ( 15.2b decay algebraically as 
n ^> oo. Moreover, these coefficients ultimatel y hav e the same sign. Thus, on the 
basis of the sufficient conditions formulated in f2071 and reviewed in Section |4] we 
can conclude that the inner jj. series in ( 14.111 ) do not converge for larger values of the 
outer index v. 

We can also insert the explicit expression (ID. 21 ) for the generalized Laguerre poly- 
nomials into infinit e seri es ( 15.21) and rearrange the order of summations, which yields 



after some algebra 112071 Eq. (3.7)]: 



Superficially, it looks as if we succeeded in constructing a power series expansion for 
the nonintegral power z^ . However, the generalized hypergeometric ser ies i Fp with 



unit argument is the limiting case z — > 1 of the so-called binomial series 111541 p. 38]: 



iFoia;z) = £ ^z'" = ^ " (-zf = i^-z)-, \z\<l. (5.9) 
A nil -^ V m J 

If we set fl = ^ — p with k e No and p e M \ No, we obtain for the iFo in ( 15.81) : 

p<0, 

iFo(^-p;l) - lim(l-z)''"* = <!o, k<p>Q, (5.10) 

z^l I 

k>p>0. 

Thus, the power series (I5.8l l for z'' is purely formal since it contains infinitely many 
series coefficients that are infinite in magnitude. 

The most general case is the Laguerre expansion (15.1b for z''e"' with p G M \ No 
and M 6 {—°°, 1/2). It is immediately obvious that z''e"*' is analytic at z = if p is 
integral, i.e., if p = m with m £ Nq, and it is nonanalytic if p is nonintegral. This 
implies that the inner ji series in (14.11b diverge if p is nonintegral, and it converges 
if p is a nonnegative integer However, p occurs on the right-hand side of ( 15.1b apart 
from the prefactor (1 — M)^"^''^'F(a + p + 1)/F(a + 1) only in the terminating hy- 
pergeometric series 2Fi{~n,a + p + l;a+l;l/(l— m)). Since the prefactor does not 
affect the convergence and existence of the resulting expansions, we have to analyze 
the asymptotics of the terminating hypergeometric series as « ^ oo. 
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Thus, with the Laguerre series coefficients 

^" -^^""^ r(a + l) 

x2/^i (-«,a + P + l;a + l;Y-7^) , peM\No, (5.11) 

the inner /i series in (14.11b diverge for sufficiently large values of the outer index v, 
and with the coefficients 

Ai"' = (l-M)""""'"'(« + l)m 

X 2^1 ( -n,« + m+l;a+l;- ), m G No (5.12) 



the inner jj. series in ( 14.1 11 1 converge. This is certainly a surprising observation which 
indicates that the difference between the almost identical terminating hypergeometric 
series 2^1 in ( 15.11b and ( 15.12b . respectively, is greater than it appears at first sight. 

Nevertheless, this puzzle can be resolved by analyzing the large n asymptotics of 
the terminating hypergeometric series in ( 15.11b and ( 15.121) . respectively. The study 
of large parameters of a Gaussian hypergeometric series is an ol d problem of sp ecial 
function theory with an extensive literature (see for example il59l4l6l[ 1 18711 and 
references therein). 

Large parameter estimates of that kind turned out be useful in the context of mul- 
ticenter integrals. In [38, Appendix], we derived and used large parameter estimates 
for some special Gaussian hypergeometric series to analyze the rate of convergence 
of certain series expansions for multicent er int egrals. 

If we apply the linear transformation 1154 p. 47] 

2Fi{a,b;c;z) - {1 -zyiFiia^c -b;c;z/{z~l)) (5.13) 

to the hypergeometric series in ( 15.121 ). we obtain: 

1 



2F1 -«,a + m + l;a + l; 



1 — M 



2Fi(-n,-m;a + l;-] , m,«eNo. (5.14) 

If n becomes large and m is fixed, the terminating hypergeometric series on the right- 
hand side can be expressed as follows: 

2Fi[-n,-m,a+l;-] = T ^ , '^\^ ^ ^. (5.15) 

The asymptotically dominant contribution on the right-hand side is the last term with 
H = m. With the help of ( 15.6b we find (-«),„ = 0(«"') as n — > 00. Since u e (-00, 1/2) 
implies m/(m— 1)g(— 1,1), we can conclude that the right-hand side of ( 15.14b decays 
exponentially as n —)► 00 because of the prefactor [u/{u~ 1)]". Thus, the inner jU series 
in ( 14.111 ) converge if the coefficients X„ are chosen according to (15.12b . 
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If we apply the linear transformation ( 15.13b to the hypergeometric series in ( 15.111 ), 
we do not obtain an expression that would be us eful for our purposes. Therefore, we 
apply instead the analytic continuation formula 111541 pp. 47 - 48] 

r(c)r(c — a — b) 

2Fi{a,b;c;z) = ^) ' \^. -^2Fi{a,b;a + b -c+l;l -z) 

r [c — aji [c — b) 

^ ria)m ^'-'^ 

X 2Fi{c — a,c — b;c ~ a — b+l;l — z), 

|arg(l -z)| < ;r, c-a-b^±m, ineNo, (5.16) 

to the hypergeometric series in ( 15.11b and obtain: 

1 



oFi -«,a + p + l;a + l;- 

1 — M 

(-P) 



2^1 a + p + l,-«;p-n + l; . (5.17) 



(a + l)„ 'V «-l 

Application of the linear transformation (15.13b then yields: 

1 

1 — M 



2F1 -n,a + p + l;a + l; 



= (1-m)"+P+' ^ ^]" 2Fi{p + l,a + p + Up-n + Uu). (5.18) 

If we set M = in ( 15.18b , the hypergeometric series on the right-hand side terminates 
after the first term. Thus, we see once more that ( 15.1b simplifies for m = to yield 
dSJl. 

With the help of ( 15.61 ) we find (p — n + 1)„, = 0(n'") as n -^ °o. Thus, the hyper- 
geometric series on the right-hand side of ( 15.18b can for arbitrary fc G No be expressed 
as follows: 

2/^i(p + l,a + p + l;p-« + l;M) 

^ (p + l).(a + p + l). ^ «^oo. (5.19) 



K=0 



{p-n + l)^K\ 



For u G (—1, 1/2), this hypergeometric series converges and we have a convergent 
asymptotic expansion as n — > 00. For u G (— °°, 1], the hypergeometric series diverges, 
but ( 15.191 ) nevertheless holds in the sense of an asymptotic expansion as n ^ 00. If 
we set ^ = in ( 15. 19b and use ( 15.71 ), we find that the Laguerre series coefficients in 
( 15.11b satisfy the leading order asymptotics 

(a) r(p + a + l) (-p)„ r(p + a + l) _p_i 
r(a + l) (a + l)„ ri-p) 

This asymptotic estimate, which does not depend on u, shows once more that z^e"' 
is not analytic at z — if p is nonintegral. The inner jj. series in ( 14.111 ) diverge for 
sufficiently large values of the outer index v. 
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Comparison of ( 15.20b with (15.7b shows that the Laguerre series coefficients X,\ 
in (15.2b . which correspond to z^e"*' with u G (— °°, 1/2), and the coefficients in (15.21 ). 
which correspond to z^ or to m = in z^e"', possess the same leading order asymp- 
totics as « — !■ oo that does not depend of u. 



6 The Transformation of One-Range to Two-Range Addition Theorems 

As discussed in the previous Sections, the legitimacy of Guseinov's rearrangement of 
a A;-dependent one-range addition theorem (13.1b . whose series expansions ( 13.101 ) for 
the angular projections ( 13.91 ) are essentially Laguerre series of the type of ( 14.4b . can 
be checked by analyzing the convergence of the transformation formula ( 14.11b . One 

only has to determine the asymptotic sign pattern and the asymptotic decay rate of 

i ex) 
the expansion coefficients X„ of the corresponding Laguerre series and employ the 

sufficient convergence criteria formulated in [202] . 

As discussed in Section |5] this approach works in a very satisfactory way in the 
case of the Laguerre series ( 15.1b for z^e"' or the equivalent one-center expansion 
(15.4b for the Slater-type function x^^(j3,r) with an in general nonintegral principal 
quantum number A^ G M \ N. 

The one-center expansions considered in Section|5]have the highly advantageous 
feature that they are comparatively simple. Therefore, we can understand the large 
index asymptotics of the expansion coefficients via the large index asymptotics of 
some Gaussian hypergeometric series, whose derivation is not straightforward but 
nevertheless not too difficult. In this way, we can explain the rate of convergence or 
divergence of the rearrangements of the Laguerre series ( 15.11 ). 

Guseinov's A;-dependent one-range addition theorems are genuine two-center prob- 
lems. Therefore, the situation is much more difficult and we are confronted with non- 
trivial technical problems. The expansion coefficients of such an addition theorem are 
according to ( 12.9b . (12.12b . and (13.1b overlap integrals, which are fairly complicated 
functions M^ ^- C whose asymptotic sign patterns and asymptotic decay rates cannot 
be determined easily. But even in this troublesome two-center case, we can arrive at 
some definite conclusions by pursuing an indirect approach based on the analysis of 
the singularities of the function which is to be expanded. 

Our starting point is a fc-dependent one-range addition theorem ( 13.11 ) for a Slater- 
type function Xni (j3 , z" ± a"') with an in general nonintegral principal quantum number 
A^ G M \ N. We assume that we succeeded in constructing the series expansions (13.10b 
for the angular projections ( 13.9b of the one-range addition theorem (13.1b . 

For the moment, let us also assume that the expansion coefficients of the series 
expansions (13.10b . which are Laguerre series of the type of ( 14.4b . have asymptotic 
sig n pat terns and asymptotic decay rates which according to the criteria formulated 
in 1I207 I1 guarantee that the transformation formula ( 14.11b produces functions C — > C 
that are analytic in a neighborhood the origin z = 0. 

Superficially, it appears that under these circumstances the legitimacy of Gu- 
seinov's rearrangement is guaranteed: A A:-dependent one-range addition theorem 
(13.1b . which is an expansion in terms of Guseinov's complete and orthonormal func- 
tions {li^ij^iY,!")} / , is transformed to a one-range addition theorem of the type of 
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( 14. Il l, which is an expansions in terms of Slater-type functions {XntiY''')} e ^i*^ 
integral principal quantum numbers n and a common scaling parameter 7 > 0. 

But such a conclusion is premature. The transformed expansion ( 14.1b in terms of 
Slater-type functions can only be one-range addition theorem, i.e., a map R^ x M^ — !■ 
C, if it converges for all r, r' e M^. This requires that the power series (14.31 ) for the 
angular projections of the expansion ( 14.21 ) for exp{Yr)XN ^iP , *" ± J"') converge for all 
r,r' e [0,00). This is a very demanding requirement, which cannot be satisfied in the 
case of exponentially decaying functions such as Slater-type functions. In the two- 
center case, it is irrelevant whether the principal quantum number of the Slater-type 
function is integral or nonintegral. 

A function /: C — >^ C is analytic at the origin z = in the sense of complex 
analysis if it has a power series in z that converges in some neighborhood of z = 0. Or 
to put it differently: Such an / is analytic at z = if its power series in z has a circle of 
convergence with a nonzero radius. As is well known, we cannot tacitly assume that 
the radius of convergence of a power series is necessarily infinite, or equivalently, we 
cannot assume that / is necessarily an entire function that is analytic for all z G C 

Therefore, we should look for features of functions analytic at z = 0, which rule 
out an infinite radius of convergence but which do not interfere with the existence 
of a Laguerre series of the type of (14.4b . My subsequent arguments are based on the 
simple, but nevertheless very consequential fact that power series and Laguerre series 
differ substantially in the way how they are affected by singularities of the function 
which is to be expanded. 

Power series converge pointwise in their circles of convergence, and in the interior 
of these circles they not only converge uniformly but they can also be used for the 
computation of higher derivatives. Since, however, the higher derivatives of a function 
ultimately become infinite in magnitude at a singularity, the radius of the circle of 
convergence is determined by the location of that singularity which is closest to the 
expansion point. 

If the function, which is to be expanded, possesses a singularity somewhere in 
the complex plane C, the radius of convergence of its power series cannot be infi- 
nite. This has immediate and undesirable consequence for integrals over the positive 
real semi-axis [0,°°) as they typically occur in the theory of Laguerre polynomials or 
in the radial parts of the three- and six-dimensional integrals in electronic structure 
theory. If the semi-infinite integration interval is not contained completely in the cir- 
cle of convergence, the term-wise integration of such a power series either leads to 
convergence to a wrong result or to divergence. 

In that respect, Laguerre series have much more convenient properties precisely 
because they in general do not converge pointwise. As discussed in more detail in 
Appendix |D] the existence of a Laguerre series of the type of ( 14.4b for some func- 
tion /: C ^ C is guaranteed as long as / belongs to the weighted Hilbert space 
^^-;,a([0,°°)) defined by dP.Sb . Of course, / must not have a non-integrable sin- 
gularity on the integration interval [0,°°), but singularities away from the positive 
real semi-axis cause no problems. Loosely speaking, we can say that a Laguerre se- 
ries simply ignores all singularities which are not located on the integration contour. 
Power series cannot do that. For them, all singularities matter even if they are far 
away from the integration contour. 
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These differences between power series and Laguerre series, respectively, can be 
illustrated by considering the following class of functions: 

F,^{z-u,e) = [z^ + u^-lzucose]"^ , z,ueC, T7,0gR. (6.1) 

For Tj = 1/2, the function Frj{z;u,9) is nothing but the explicit expression for the 
difference |jr — j?! = [x^+y^ — 2xycos9Y'^ of two vectors j:,>' G K^ in spherical polar 
coordinates in disguise, and for Tj = — 1 /2, it corresponds to the Coulomb or Newton 
potential l/jjc — j'j. 

As long as F-q{z\u,Q) does not have a non-integrable singularity on the posi- 
tive real semi-axis, it belongs to the weighted Hilbert space L^--.a ([0, °°)) defined in 
(ID. 8b . Therefore, Fri{z;u,9) possesses a Laguerre series of the type of (14.4b . although 
I was not able to find a closed form expression for the coefficients of this expansion. 

If TJ is a nonnegative integer, rj ~ n with n € No, Fr^{z;u,9) is a polynomial in 
z and therefore an analytic function for all z e C But if Tj e R \ No, Fr^{z;u,9) has 
singularities at 

zi,2 = |cos0±y^[cos0]2-l|M. (6.2) 

Accordingly, the radius of convergence of its power series about z = is equal to \u\. 
The function Frj{z;u,9) essentially corresponds to the well known generating 
function 

[l-2xf + f2]-^ = £c„^(x)f«, |f|<l, A^O, (6.3) 

n=0 



of the Gegenbauer polynomials 0154 p. 222]. Thus, F-q{z',u,9) possesses - depending 
on the relative magnitudes of |z| and \u\ - two complementary power series expan- 
sions. For \z/u\ < 1, it possesses a convergent power series in z/u, 

Fri (z; M, 0) = M^'' £ C-l (cos 0) {z/uf , (6.4) 

«=o 

and for \z/u\ > 1, it possesses a convergent inverse power series in z/u, 

F^ (z; «, 0) = z^" £ C-" (cos 0) («/z)" , (6.5) 

«=o 

which can also be interpreted as a convergent power series in u/z- 

Let us now assume that the singularities ( 16.21 ) of F-q{z;u,9) with 77 G R \ N do 
not lie on the integration interval [0,°°) and that we succeeded in constructing its 
Laguerre series of the type of (14.4b . possibly by purely numerical means. Then, the 
application of the transformation formula (14.11b to this Laguerre series produce s - if 
necessary with the help of nonlinear sequence transformations as described in 12071 
Section 6] - the power series (16.4b . which converges for \z/u \ < 1 and which diverges 
for \u/z\ < 1. 

The small z series (16.4b accomplishes at least for \z/u | < 1 a separation of the 
variables z and u and therefore resembles a two-range addition theorem. For \z/u\ > 1 , 
the large z series ( 16.5b also accomplishes this separation, but I see no obvious way of 
computing the large z series (16.5b from the Laguerre series for F,j (z; m, 0). Instead, we 
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would have to construct a Laguerre series in u - this is possible since z and u play a 
symmetrical role in F,j (z; m, 0) - from which we can compute the large z series ( I6.5l l 
which is also a small u series. 

These considerations show that if we want to represent Frj{z;u,Q) by power se- 
ries in z, a two-range scenario cannot be avoided. It does not matter if we start from 
a Laguerre series which provides a unique representation of F-q (z; m, 0) that is com- 
putationally useful in integrals over the whole real semi-axis. The singularities ( 16.21 ) 
of Fri{z',u,9) with Tj e M\N imply that there can be no power series in z which 
converges for all z G C 

In the context of addition theorems, it may be of interest that the generating func- 
tion ( 16.31 ) of the Gegenbauer polynomials can be used for the construction of two- 
range addition theorems in a relatively straightforward way. If we set A = — v/2, 
t — r^/ry, and x = cos0 in the generating function (16.3b . we obtain the following 
Gegenbauer expansion for the general power function 



|r<±r>r -r^£(Tl)"Cr/'(cos0)(r</r>)«, veM. (6.6) 

n=0 



The two-range form of this Gegenbauer expansion is a direct consequence of the 
convergence condition \t\ < 1 in (16.31) which translates to the convergence condition 
r</r> < 1. 

We can easily construct an addition theorem from the Gegenbauer expansion 
( 16.6b : We only have to replace the Gegenbauer polynomials by Legendre polyno- 
mials. However, the practical realization of this obvio us idea had apparently not been 
so easy. As discussed by Steinborn and Filter |182, pp. 269 - 270], many authors 
had quite a few problems with the determination of explicit expressions for the co- 
efficients of the expansion of Gegenbauer polynomials in terms of Legendre poly- 
nomials. Also Steinborn and Filter constructed very messy expressions for these co- 
efficients which are restricted to certain superscripts of the Gegenbauer polynomial 



1 1811 Section 31. 

This is somewhat strange because already at that time a much more convenient 
expression for these expansion coefficients had been a vailab le in the mathematical 
literature. In Exercise 4 on p. 284 of Rainville's book 116511 . one finds the follow- 
ing relationship, where [m/2j denotes the integral part of m/2 (compare ll218l Eq. 
(5.2)]): 

Ci; (x) = ^'"f ^^Vdf'V^^^' (2— 4. + 1 ) P,.-2s{x) . (6.7) 

.v=0 [J/^)m-sS. 

This result can be proved via the explicit expression [36, Eq. 7.313.7 on p. 836] for 
the integral j\{l -x)"{l -xy-^/^Ci!r{x)C^{x)dx. One only has to set v = 1/2 and 
perform the limit a ^ 0, which requires, however, some algebraic trickery. 
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If we now insert ( I6.7l i into (16.6b and rearrange the order of summations, we obtain 
after some algebra the following expansion in terms of Legendre polynomials: 

., (-V/2), 



(3/2)^ 



2Fi(£-v/2,-[v+l]/2;^ + 3/2;[r</r>]2). (6.8) 



If V is an even integer, v = 2n with n G No, |r< ± r> |*' is a polynomial in both r< 
and r< and therefore analytic. The infinite i series in (16.81 ) terminates because of 
the Pochhammer symbol (— v/2)^ = (—«)<', implying £ <n. Similarly, the Gaussian 
hypergeometric series 2F1 in ( 16.8b terminates since £ — v/2 ^£ — nis either a negative 
integer or zero. 

In (16.8b . we only have to replace the Legendre polynomials by spherical harmon- 
ics via the so-called spherical harmonic addition theorem ( 12.4b to obtain an expa nsion 
in terms of spherical harmonics, which had originally been derived by Sack il69l 
Eq. (19)] by solving a partial differential equation and which converges as long as 
r</r> < 1 holds; 

i=Q m=-i 

X ^-^2^i(^-v/2,-[v+l]/2;£ + 3/2;[r</r>]2). (6.9) 
This two-range addition theorem simplifies considerably and also assumes a one- 



range form if v is a positive even integer, v = 2n with n G N (see above or also 11202 , 
pp. 1258 - 1259]). But for arbitrary v G E, ( 16.91 ) is a two-range addition theorem. 
This is a direct consequence of its derivation via the generating function ( 16.3b of the 
Gegenbauer polynomials, which is a power series in t with a nonzero, but finite radius 
of convergence. 

This observation suggests the following interpretation: A two-range addition the- 
orem for a function /(r< ± r>) having its only singularity at r< ± r> =0 corre- 
sponds to a possibly rearranged power series in r< with a finite radius of conver- 
gence, and this radius is determined by the condition r< < r>. This interpretation is 
also confirmed by the differential operator ( 12.71 ) which had been the central tool for 
the derivation of two-range addition theorems in 1 199, 201] or in i202t, Section 7]. 

This interpretation also applies to the addition theorem (16. 91) . For arbitrary v 7^ 
0,2,4,..., the general power function |r< ±r>|*' has a singularity for r^zLr^ — 0, 
and this singularity determi nes the radius o f convergence of its power series in r< and 



enforces a two-range form 111 99112011120211 . 



These conclusions about the nature of the addition theorem for the general power 
function |r< ± r>|^ help us to understand some essential features of addition theo- 
rems for Slater-type functions that converge pointwise. For the sake of simplicity, let 
us first consider the so-called Is function 

exp(-)3|r-r'|) = exp(-j3[r2 + r'2-2rr'cos0]i/2) (6.10) 
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and let us also assume r' > 0. If Guseinov's rearrangements are legitimate, then 
exp(— )3 |r — /•'!) must possess an expansion of the type of ( 14.11 ) in terms of Slater- 
type functions {x'"£{Y, ^)} t ' which converges for the whole argument set M? x K.^. 
Obviously, this is equivalent to requiring that exp(7r) exp(— j3 \r — r'|) possesses a 
power series about r = 0, which converges for all r e [0,oo). Since, however, exp(7r) 
is an entire function, whose power series in r converges in the whole complex plane 
C, we can ignore it for the moment. It is sufficient to analyze whether the Is function 
( I6.IOI 1 possesses a convergent power series in r and for which values of r and r' this 
series converges. 

If we expand the exponential on the right-hand side of (16.101 ). we obtain: 

exp(-j3|r-/-'|) = £ ^^^^[r^ + r'^-lrr'cosO]"/^. (6.11) 

If the index Kis even, K' = 2^ with A; e No, then [r^ + r'^ — 2rr' cos 9]'' is a polynomial 
in r which is obviously analytic for all r G [0,oo). But if K is odd, K ~2k+l with 
k e No, we are in trouble. We can use the generating function (16.3b of the Gegenbauer 
polynomials to obtain a power series expansion in r/r' 



..2 , /2 



r 



-2rr'cos0]'^+'/2 ^ /2'^+i [1 + (^//)2 _2(r/r')cos 0]'+'/' 



j2k+l v^ ^-k-l/2 



l^Cj'-'''{cos9)ir/r'y. (6.12) 

;=o 

Unfortunately, the Gegenbauer expansion ( 16.12b converges only if r/r' < 1 holds. 
Therefore, the expansion obtained in this way corresponds to the case \r\ < \r'\ of the 
two-range addition theorem for |r ± r^^''^ which is a special case of (16. 9b . 

This has some far-reaching consequences for the analyticity of the Is function 
exp(— )3|r — r'l) at r = 0. The odd powers in the power series expansion (16.111 ) 
are singular for r — r' — 0. This fact makes it impossible to construct for either 
exp(— /3|r— r'l) or exp(7r)exp(— /3|r— r'|) with r' > a power series in r that con- 
verges for all r £ [0,°°). This is only possible in the one-center case r' = 0. Conse- 
quently, a one-range addition theorem of the type of (14.11 ). which converges point- 
wise for all r G M? or equivalently for all r G [0,°°), cannot exist for the Is function 
ifr'^0. 

An explicit expression for two-range addition theorem for the \s function (16.101 ) 
can be derived easily via the following Gegenbauer-type addition theorem for the 

modified Bessel function w^^K^ [yw) with w = [p^ + r^ — 2rp cos 0] , < p < r, 
and V G C\No L154^ pp. 106 - 107]: 

w-^Kyirw) = 2^7-^r(v)(rp)-^ f; C;(cos0)/v+„(7p)^v+„(7'-)- (6.13) 

Here, Iv+n{YP) an d Ky +niyr) are modified Bessel function of the first and second 



kind, respectively 11154 p. 66]. 

The modified Bessel function w^^Ky (pv) in ( 16.131 ) is essentially a reduced Bessel 
function ky (yw) defined by ( 12.16b . On the basis of ( 12.41 ). ( 16.7b . and ( 16.13b . the follow- 
ing two-range addition theorem for reduced Bessel functions with half-integral orders 
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can be derived in a fairly straightforward way r ill82l Eq. (3.4)] or, as an improved ver- 
sion 12181 Eq. (5.5)]): 



^„_,/2(i3|r<±r>|) 



i^^(j3r<)"-i/2()3,>)''-i/2 



(2n-l)!! 
I I (Tl)'^[i'r(r<A<)]*Ff('->A>) 



e=o m=- 

" (-«)v(l/2-«),^ 
,^0 v!(3/2),+v ^ 



-2v- 



1/2) 



X k+2v-n+l/2{Pf<)Kt 



i+2v-n+l 



/2(j3r>). 



(6.14) 



This addition theorem was quite consequential for my later scientific interests. In my 
diploma thesis 11921 . which was published in condensed form in 11831 . I used this 
addition theorem for the evaluation of simple multicenter integrals of reduced Bessel 
functions. 

If we set m = in ( 12.171 ). we obtain ^1/2(2) = e^^". Thus, we only have to set « = 1 
in (16.14b to obtain a two-range addition theorem for the Is function. At first sight, 
neither the addition theorem ( I6.14l i nor its special case with n ~ I looks like a power 
series expansion in r<. However, the m odified Bessel functions Ie+2v-n+i/2{P''<) is 
defined by a power series in /3r< 11154 p. 66], which shows that the addition theorem 
(16.14b is nothing but an infinite multitude of ^-dependent rearranged power series 
expansions in r<. 

The same conclusions hold for the following two-range addition theorem of a B 
function [201, Eq. (5.11)] which can be viewed to be an anisotropic generalization of 
the addition theorem (16.14b for a reduced Bessel function: 



^:'(r<) 



\ ^2 (Rr \n+l-li-q-ll2j 

——[lir<) ^ in+e+ei-q+l/2 



^oin + i-q) 



iPr^) 



V ^ ^ (£2111 -\- mi\iim\\im) 



s=0 



A £2 



B"'+'"' (B 



>)■ 



(6.15) 



Other two-range addition theorems of B functions are discussed in 112181 Sections IV 
andV]. 



In (16.15b . {£2m + mi \£\m\ \£m) is a so-called Gaunt coefficient 113 If which corre- 
sponds to the integral of the product of three spherical harmonics over the surface of 
the u nit sphere in R^. The selection rules of this Gaunt coefficient (see for example 
I 21I Section 3] or |202, Appendix C]) imply that A£2 = (^ + ^1 -^2)/2 in (l6J3T l is 
always either zero or a positive integer. The symbol Y/ in ( 16.15b indicates that the 
summation proceeds in steps of two. 
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Since a Slater-type function with an integral principal quantum number can ac- 
cording to ( 12.201 ) be expressed as a finite sum of B functions, we can conclude that 
the two-range addition theorem of a Slater-type function obtained by forming linear 
combinations of (16.15b is nothing but an infinite multitude of rearranged power series 
expansions in r< . 

These considerations apply also to Slater-type functions ;i^^^(/3,r< ±r>) with 
nonintegral principal quantum numbers A^ or to exp(7r);|fj^^ (j3 , r< ± r> ) . These func- 
tions are obviously singular for r^±r^ = 0, which implies that their power series 
expansions about r< = can only converge for r< < r>. Consequently, an expansion 
of the type of (14. Il l in terms of Slater-type functions {x^\{7 >'')}». Cm with integral 
principal quantum numbers « and a common scaling parameter 7 > 0, that converges 
pointwise for all r, r' G M^, cannot exist. 

These considerations can be generalized further: Let us assume that /(r ± r') with 
/ > is singular for r ± r' = but analytic elsewhere. Accordingly, a power series 
expansion for either /(r ± r') or for e'^'' f{r ± r') about r = can only converge for 
r < r'. This rules out the existence of an expansion of f{r±r') in terms of Slater- 
type functions X'^eiYi'')} with integral principal quantum numbers n and a common 
scaling parameter 7 > 0, that converges pointwise for all r, r' G M^. 

As shown in 11 1991120 11 or in [202, Section 7], pointwise convergent addition the- 



orems are nothing but rearranged Taylor expansions. Thus, the assumed singularity of 
such an /(r ± r') for r±r' = implies that a pointwise convergent addition theorem 
must have a two-range form. 

Slater-type functions as well as all the other commonly used exponentially de- 
caying basis functions have a singularity at the origin. Consequently, their pointwise 
convergent addition theorems must have a two-range form. In contrast, the Gaussian 
function exp(— j3|r±r'p) is analytic for all r,r' G R ^. Co nsequently, it possesses a 
one-range addition theorem that converges pointwise 11 144 , Eq. (8)]. 

It makes no difference if we start from a Laguerre series for a given function. If we 
apply the transformation formula ( 14.111 ) to the series expansions ( 13.101 ) for the angular 
projections ( 13. 9t of Guseinov's one-range addition theorem ( 13.1b for x^^(j3, r±r'), 
we obtain power series expansion in r for the angular projections which converge or 
diverge, depending on the relative magnitude of r and r'. 

Let us now assume that /(r ± r') is a function which is singular for r ± r' = 
and that we know its two-range addition theorem. If we accept the premise that its 
two-range addition theorem is nothing but a multitude of rearranged power series 
expansions in r that converge for r/r' < 1, then the uniqueness of a power series in 
the interior of its circle of convergence implies that Guseinov's rearrangements of 
one-range addition theorems cannot produce anything new. We obtain a two-range 
addition theorem for exp(7r)/(r ± r') by forming the Cauchy product of the power 
series for exp(7r) with the possibly rearranged ^-dependent power series in r that 
occur in the two-range addition theorem for /(r ± r'). In the final step, we only have 
to multiply the resulting addition theorem for exp(7r)/(r ± r') by exp(— 7r) or its 
power series to arrive at the possibly rearranged addition theorem for /(/• ± /) from 
which we started. 
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The discussion of this Section may create the false impression that singularities 
of basis functions for electronic structure calculations are something very negative. 
This is certainly true in the case of analytical manipulations since virtually all manip- 
ulations becomes more difficult in the presence of singularities. However, Kato 1114111 
had shown that the singularities of an atomic or molecular Hamiltonian translate to 
corresponding singularities of the eigenfunctions, commonly called cusps. But a cusp 
is just another word for a singularity. Thus, the ability of basis functions to reproduce 
the singularities of exact wave functions is of considerable importance for the rate of 
convergence of an electronic structure calculation. 

Gaussian functions do not have singularities like exponentially decaying func- 
tions. In my opinion, this is the main reason why their multicenter integrals can be 
evaluated much more easily than the corresponding integrals of exponentially decay- 
ing functions. At the same time, the absence of singularities is also a major drawback 
of Gaussian functions. Many Gaussian functions are needed to approximate functions 
possessing singularities with sufficient accuracy. 



7 Numerical Implications of Truncated Expansions 

The analysis of Sections[T]and|5]shows that the one-center expansion (11.2b of a Slater- 
type function XNiiP,!") in terms of Slater-type functions {^^^(7, >■)} _, , , with in- 
tegral principal quantum numbers « e N and a common scaling parameter 7 > does 
not exist if the principal quantum number A^ is nonintegral. The leading coefficients 
of (11.2b are zero, and the higher coefficients are infinite in magnitude. 

The nonexistance of (11.2b is a direct consequence of the fact that the radial part 
of XNiiP,!") with A^ e M\N is not analytic at r = 0. As shown in Section |5j this 
nonexistance can also be shown by applying the transformation formula ( 14.11b to the 
radial part of a fc-dependent one-center expansion ( 15.4b for Jj^^(j3,r). In the case 
of nonintegral quantum numbers, the transformation formula (14.111) produces power 
series coefficients that are either zero or infinite in magnitude. 

In Section |6] it was shown that the rearrangement of a A:-dependent one-range 
addition theorem (13.1b . which is a much more complicated expansion than its one- 
center limit (15.4b . also does not produce the desired result. Both ;|f^^(j3,r ± r') and 
exp(7r);i^j^^(/3,r± r') with r' > are analytic at r = 0. Consequently, the appli- 
cation of the transformation formula ( 14.11b to the series expansions of the angular 
projections of exp{Yr)Xi!j lifi i^^f^) produces mathematically meaningful power se- 
ries expansions in r. Unfortunately, these power series expansions have a finite radius 
of convergence and they converge only for r < /. This follows at once from the fact 
that both XNiiP j *■ =•= /"') and exp(7r);^^^()3 , r ± r') are singular for r±r' — 0. Ac- 
cordingly, Guseinov's rearrangements transform a one-range addition theorem for 
xltfi{li,r±r') to a two-range addition theorem. This is certainly not what Guseinov 
had tried to achieve. 

But these observations do not tell the whole truth, and in particular they do not im- 
ply that Guseinov's approach is necessarily doomed. There is overwhelming evidence 
that Guseinov never transformed a complete one-range addition theorem of the type 
of ( 13.11 ) containing an infinite number of terms. Instead, Guseinov only transformed 
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truncations of either addition theorems or their one-center limits, which all contain a 
finite number of terms only. Thus, Guseinov exclusively did his rearrangements with 
an analog of ( I4.8l l, which transforms a truncated Laguerre series of the type of ( 14.6b to 
a polynomial of the type of ( 14. 9t . To the best of my knowledge, Guseinov never used 
the complete transformation formula ( 14.111 1 and only claimed - probably on the basis 
of an insufficient amount of numerical evidence - that his transformed truncations 
remain meaningful in the limit of infinite expansion lengths, but he never provided 
convincing evidence supporting his claim. 

The radial parts of the angular projections of Guseinov's truncations are apart 
from a common exponential exp(— yr) finite linear combinations of generalized La- 
guerre polynomials in 2yr. Thus, for finite truncation orders JV , Guseinov's rear- 
rangements are legitimate and produce polynomials in r. 

In actual calculations, we always have to truncate infinite series expansions after 
a finite number of terms unless we are fortunate enough to find a way of expressing 
a series in closed form. Therefore, a skeptical reader might argue that Guseinov's 
approach is completely satisfactory from a practical point of view, and that my insis- 
tence on the convergence and existence of infinite series expansions, which in actual 
calculations have to be truncated and thus do not really occur in practice, is nothing 
but a mathematical over-sophistication. 

However, Guseinov's rearranged addition theorems do not only suffer from the 
fact that in limit of infinite expansions lengths they either do not exist or that they 
lose their one-range nature. These rearrangements cause also other problems which 
become particularly evident if one tries to achieve a (very) high accuracy by including 
a large number of terms in Guseinov's truncations. 

One serious problem was already discussed in Section [3] As explained there, I 
have grave doubts that a Guseinov function ^f^™ (7, t) can be expressed in a numeri- 
cally stable way as a finite sum of Slater-type functions {z™^(7>'')}„/=^ + 1 ^i^ (13.21 ) 
if the index n is large. This applies also to the substitution of overlap integrals in- 
volving Guseinov functions by finite sums of overlap integrals involving Slater-type 
functions via (13.3b . Both ( 13.2b and ( 13.3b are based on the explicit expression ( ID. 21 ) 
of the generalized Laguerre polynomial L„ (z), which tends to become numerically 
unstable if its index n becomes large. The reason is that the coefficients of orthogonal 
polynomials have strictly alternating signs. 

The substitution of a function set, which is complete and orthonormal in a given 
Hilbert space, by a function set, which is only complete, but not orthogonal, also 
causes some nontrivial problems. In Section[3]it was emphasized that orthogonal ex- 
pansions tend to be computationally well behaved because Parseval's equality ( 13.61 ) 
guarantees that their series coefficients are bounded in magnitude and vanish for large 
indices. In contrast, the series coefficients of nonorthogonal expansions are not nec- 
essarily bounded in magnitude and do not necessarily vanish with increasing index. 

These complications with unbounded coefficients can be demonstrated convinc- 
ingly by considering a truncation of the comparatively simple one-center expansion 
( 15.4b . This ansatz corresponds to the following approximation of a Slater-type func- 
tion XfjiiP ' *") ^i'^h an in general nonintegral principal quantum number A^ G K \ N 
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by a finite sum of Guseinov functions: 



XffAP^r) 



[li + Y]N+L+k+2 (2L + k + 2 



v=0 



V! 



A.fv'tL+i.Llr,'-) 



X2Fi(-V,N + L + k + 2;2L + k + 3--^] . (7.1) 

V p + yJ 

If XNiiP^'") belongs to the weighted Hilbert space L\{R^) defined by ( 12.151 ). this 
•vK-dependent approximation is mathematically meaningful, and we can be sure that 
it converges in the mean with respect to the norm of L\(K^) as JK -^ o°. 

We only have to cancel the spherical harmonics and remove the exponential on the 
right-hand side of (17.11 ) to see that the finite sum on the right-hand side corresponds 
to a truncated Laguerre series of the type of ( 14.61 1. If we now apply the transformation 
formula (14.81) to this truncated Laguerre series, we ultimately obtain an approximation 
of the Slater-type function Xnl (J3 , *■) with A^ G M \ N as a finite linear combination of 
Slater-type functions X^jjyJ-r) with integral principal quantum numbers n. 

For finite values of the truncation order jV , all coefficients of this approximation 
to x^ ijfi^r) are well defined and finite. Unfortunately, this does not imply that this 
expression remains well behaved in the limit jV —> oo. For large values of JV , the 
leading coefficients of the resulting expression approach zero and the higher coeffi- 
cients diverge in magnitude. 

This can also be demonstrated by inserting the expansion coefficients of the La- 
guerre series ( 15.21 ) for z^e"^' with nonintegral p G M \ No into the \i series in ( I4.11l i. If 
we replace the exact expressions for the Laguerre series coefficients of the later terms 
of the jU series with indices /x > M by their leading order asymptotic approximations 
( 15.20b . we obtain the following leading order asymptotic approximation: 

The infinite series on the right-hand side is nothing but the tail of the Dirichlet series 
^{^) = E r-n { n + 1)^* for the Riemann zeta function with 5 = p — v + 1 (see for 
example 01541 p. 21]). As is well known, the Dirichlet series for C,{s) converges for 
5R(.s) > 1 and diverges for '^{s) < 1. Thus, the series (I7.2l i converges for V < p, and 
it diverges for v > p . 

If 3i{s) is only slightly larger than 1, the convergence of the Dirichlet series 
Lr=o('' + 1) * can become prohibitively slow (the horrifying example of the Dirich- 
let series with s = 1.01 is discussed in 12141 p. 194]). The slowest convergence on 
the right-hand side of ( 17.21 ) occurs for the largest value of v satisfying V < p, i.e., for 
V = [pj, where [pj is the integral part of p. In this case, convergence can become 
so slow that it is practically impossible to evaluate the series on the right-hand of 
( I7.2l i with sufficient accuracy by adding up its terms successively. Instead, one has 
to employ suitable convergence acceleration techniques as for example the Euler- 
Maclaurin formula (see for example 112141 Section 2] and references therein). In this 
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context, it may be of interest that the Euler-Maclaurin formula for the Riemann zeta 
function and other asymptotic approximations to truncation errors of series represen- 
tations for special functions can also be derived by solving systems of linear equations 
i203ll . 

A skeptical reader might argue that my digression on the convergence properties 
of the Dirichlet series for the Riemann zeta is of no interest in the context of molec- 
ular electronic structure. However, I personally became interested in the Riemann 
zeta function because of certain infinite series expansions occurring in expressions 
for molecular integrals of exponentially decaying functions involving the Coulomb 
potential. The convergence properties of these expansions closely resemble that of 
the Dirichlet series for l^{s) with 3i{s) not much larger than 1 (see for example 1381 
Tables I, II, V, VI, and VII] or |184, Table 1]). 

The slow convergence of the series on the right-side of (17.2b for v only slightly 
smaller than p makes it (very) hard or even practically impossible to observe by 
purely numerical means that the leading terms of the formal power series for z^s"'' 
with nonintegral p G M \ No vanish. 

Moreover, the series on the right-hand side of ( 17.2b diverges for v > p, but it 
diverges quite slowly if v is only slightly larger than p . Quite a few terms are needed 
to observe the divergence of such an infinite series in the case of small or moderately 
large values of v. There is also the additional complication that it is by no means 
easy to establish unambiguously the divergence of a series of the type of ( 17.2b if only 



a finite number of terms are available 1115811 . 



A further complication occurs if we do not rearrange an infinite Laguerre series 
of the type of ( 14.4b via the complete transformation formula ( 14.11b . which uses all 
Laguerre series coefficients A„ with n G No, but via ( 14.8b which transforms the 
partial sum /m (z) of the Laguerre series ( 14.4b to a polynomial of degree Minz- Since 
the transformation formula (14.8b only uses the coefficients Xq , Aj , . . . , A^ and 
since M is in practice only moderately large, neither the vanishing of power series 
coefficients with v < p nor the divergence of the coefficients with v > p can be 
observed easily (everything remains finite). Consequently, I would not be surprised if 
Guseinov and Mamedov, who certainly had don e test calculations, simply overlooked 



the nonexistence of their one-center expansion 111221 Eq. (4)] in the Umit ^ — > °°. 

It would have been interesting if Guseinov and Mamedov had applied sequence 
transformations, because this could have helped them to see the vanishing or the 
divergence of their coefficients more clearly. My suggestion may sound paradoxical 
because sequence transformations are normally used to accelerate convergence or to 
associate a finite value to a divergent sequence or series. It was, however, shown in 
recent articles by Beckermann, Kalyagin, Matos, and Wielonsky (81], Beckermann, 
Matos, and Wielonsky |9], Brezinski Ill6ll . Brezinski and Redivo Zaglia |17], and 
Guilpin, Gacougnolle, and Simon [39'| that sequence transformations can also be used 
to determine the location of discontinuities of functions more precisely or to show 
them more clearly. By a slight abuse of language, such an application of sequence 
transformations could be called acceleration of divergence. 

Guseinov's and Mamedov's inability of observing any problems with their rear- 
ranged one-center expansions highlights once more the dangers of relying entirely on 
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numerical test calculations without trying to understand the subtleties of the underly- 
ing mathematics. 

One should also take into account that the apparent convergence of the sum of 
the leading terms of an infinite series to the correct limit does not prove the existence 
of this series, let alone its converges to the correct limit. As discussed in Appendix 
IE] the phenomenon of semiconvergence - initial apparent convergence of the leading 
terms of an infinite series followed by divergence if more terms are included - is well 
established in the literature. 

The examples in Appendix|E]on semiconvergence and related phenomena should 
suffice to convince even a skeptic that misinterpretations of seemingly convincing 
numerical evidence can happen easily. It also happened to me. In Ill96tl I misinter- 
preted my su mmat ion results obtained by applying Wynn's epsil on algorithm 1219], 
the d variant il95i Eq. (7.3-9)] of Levin's transformation 115111 . and the d variant 
of the so-called -Y transformation fl95, Eq. (8.4-4)] to the partial sum of the first 
22 perturbation series coefficients of the factorially divergent Rayleigh-Schrodinger 
perturbation series for the ground state energy of the quartic anharmonic oscillator 
All calculations were done in FORTRAN 77 on a Cyber 1 80-9 95 E with a precision 



of approximately 29 decimal digits (for more details, see 112051 pp. 7 - 8]). 

Unfortunately, my conclusion in lll96ll that Levin's transformation produces a 
convergent result - although perfectly plausible at that time - was premature and 
based on th e inc omplete evidence provided by the first 22 perturbation series coef- 
ficients. In 121311 . we repeated my previous calculations using now 200 perturbation 
series coefficients calculated exactly with the help of Maple's rational arithmetic, and 
we did the summation calculations in Maple with a precision of up to 1000 decimal 
digits and transformation orders as high as ^ = 199. These calculations showed un- 
ambiguously that Levin's transformation failed to produce convergent results in the 
case of higher transformation orders, and that this failure could not be attributed to 
numerical instabilities. The divergence of Levin's transformation was also confirmed 
in 1 197, Table 2]. A similar divergence of Levin's transformation was later observed 



by Cizek, Zamastil, and Skala 11191 p. 965] in the case of the hydrogen atom in an 
external magnetic field. 

Of course, my misinterpretation in fl96T or similar problems of other authors do 
not rule out the possibility that carefully conducted purely numerical investigations 
can provide valuable theoretical ins ight. In llioll we formulated with the help of some 
numerical techniques developed in 120011 the conjecture that the factorially divergent 
perturbation expansion of a certain non-Hermitian ,^,5^-symmetric anharmonic os- 
cillator is a Stieltjes series. Recently, our conjecture, whose correctness implies the 
Fade summability of this perturbation expansion, was proved rigorously by Grecchi, 
Maioli, and Martinez |37]. 

Let us now assume that we rearrange a truncation (17.1b of the one-center ex- 
pansion ( 15.4b for xiviiP I *") ^i'^h nonintegral principal quantum number N eM.\N. 
If the truncation order ^ in ( 17.11 ) is small, its increase will certainly improve the 
accuracy of the rearranged truncation. However, for sufficiently large values of ^, 
the behavior of the rearranged truncations changes. The leading order asymptotic 
approximation ( 17.2b implies that the accuracy ultimately deteriorates with increas- 
ing ,yK, and for -yV -^ °a, the rearranged truncations ultimately become mathemati- 
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cally meaningless. Thus, rearranged truncations of the one-center expansion ( I5.4l i are 
semiconvergent with respect to a variation of JV for nonintegral principal quantum 
numbersA^eR\N. 

The semiconvergence of the rearranged truncations implies that they can be used 
for computational purposes at least for sufficiently small truncation orders jV . But 
obviously, one should be careful. It is necessary to investigate for which values of 
jY the intrinsic pathologies of the rearranged expansions become intolerable. The 
leading order asymptotic approximation ( 17.21 ) indicates that even fairly large values of 
jY should produce acceptable results, but additional and in particular more detailed 
numerical investigations, which also try to estimate the detrimental effect of possible 
numerical instabilities, would certainly be desirable. 

One can look at the rearrangements of truncations ( 17.11 ) also from a different per- 
spective. As discussed in AppendixICl finite approximations to a function / e ^ of 
the type of dC.ll ) in terms of a function set, that is complete but nonorthogonal in 
some Hilbert space M', can be constructed by minimizing the mean square devia- 
tion ( IC.2b . although we cannot tacitly assume that these finite approximations can be 
extended to infinite expansions of the type of (IC.3I) . 

Thus, is certainly legitimate to approximate a Slater-type function XiviiP i *") ^^'^^ 
a nonintegral principal quantum number A^ S M \ M by a finite linear combination of 
Slater-type functions X^iXy, r) with integral principal quantum numbers n. However, 
a determination of the coefficients of these approximations via a purely numerical 
minimization of the mean square deviation does not look like a promising computa- 
tional strategy in the context of addition theorems or multicenter integrals. In such a 
case, Guseinov's approach is most likely the better alternative. 

The truncations ( 17.1b of the one-center expansion ( 15.4b are simple enough to per- 
mit a detailed mathematical analysis of Guseinov's rearrangements. But we would of 
course be much more interested in understanding the subtleties of the rearrangements 
of the truncations 



II I 



Z*^,(i3,r±r') « I I I ,x;:;^i'(7,^,±r'),f^;",(7,r) (7.3) 



of a A^-dependent addition theorem ( 13.1b . The overlap integral i^„j,'„^ (/jj^ji'"') in 
(ITJT i is defined by by dTTB . 

Since the A:-dependent addition theorems (13.1b are genuine two-center problems, 
it is not at all easy to do a rigorous analysis of the mathematical properties of re- 
arrangements of the truncations (17.3b . Let me emphasize once more that there is a 
fundamental difference between the one-center case analyzed in Section |5] and the 
two-center case analyzed in Section |6] In the limit of infinite truncation orders jY , 
Guseinov's rearrangements of one-center expansions for Slater-type functions with 
nonintegral principal quantum numbers produce mathematically meaningless expan- 
sions, and for finite values of ,JY these rearrangements are semiconvergent. 

In the two-center case, Guseinov's rearrangements of one-range addition the- 
orems produce mathematically meaningful expansions, but they only converge for 
r < r'. Consequently, these expansions correspond to the small r parts of two-range 
addition theorems. This is certainly not what Guseinov had tried to achieve. 
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Apart from being an undesirable result, the two-center nature of rearrangements 
of the truncations (17.3b in the limit of infinite truncation orders .jV is also a possible 
source of problems. For finite values of jV , the truncations ( I7.3l l are obviously one- 
range addition theorems: The vector r occurs exclusively in the Guseinov functions 
j,f^'" (7, r), and r' occurs exclusively in the overlap integral f^,^\'„^ (7ij3,±''')- This 
applies also to their rearrangements. The vectors r and r' are still separated even if 
the Guseinov functions are replaced by Slater-type functions with integral principle 
quantum numbers via ( I3.2l i and the overlap integrals involving Guseinov functions 
by overlap integrals of Slater-type functions via (I3.3l l. 

But in the limit of infinite truncation orders J/ , a change resembling a phase 
transition takes place: The resulting rearranged expansions lose their advantageous 
one-range nature and converge to two-range addition theorems. As is well known, 
the decay rate of the coefficients 7„ of a power series /(z) — L"=o 7nZ" as « — > 00 
determines the convergence type of this series: If 7, decays factorially as « ^^ 00, f[z) 
is entire, and if 7„ decays only exponentially, the radius of convergence of its power 
series is finite. Thus, the expansion coefficients of the Slater-type functions x'ndYj'') 
in Guseinov's rearranged addition theorems apparently decay at most exponentially 
since these expansions only converge for r < /. 

Otherwise, very little can be said about the ^-dependence of Guseinov's rear- 
rangements of the truncations ( 17.31 ). The problem is that in Section|6]the mathematical 
properties of Guseinov's rearranged addition theorems were analyzed via the singu- 
larity of a Slater-type function xivLiP ,i'±r') for r±r' ^ 0, but not via the trans- 
formation formula ( 14.111 ). We would need the large index asymptotics of the overlap 
integrals ^X^ y^^ {Y,P, i'"') in ( 13.11 ). On the basis of our current level of understand- 
ing, such an asymptotic analysis seems to be out of reach. We do not have anything 
resembling a two-center analog of the leading order asymptotic approximation ( 17.21 ) 
which turned out to be so very useful in the one-center case. 

Guseinov's rearrangements of one-range addition theorems can be used to con- 
struct approximations to a Slater-type function XNiiP )'"='=''') with both integral and 
nonintegral principal quantum numbers in terms of a finite number of Slater-type 
functions Z™£(7, r) with integral principal quantum numbers n. As in the one-center 
case, such a finite approximation of the type of (IC.ll ) can - at least in principle - be 
constructed by minimizing the mean square deviation ( IC.2b . although such a numeri- 
cal determination of the expansion coefficients is most likely not a very good idea in 
the context of multicenter integrals. 

But there remains a principal problem. If we use a two-range addition theorem in 
a multicenter integral and do not take into account its two-range nature by splitting up 
the integration interval of the resulting radial integration, we may well end up either 
with convergence to the wrong Umit or even with a divergent series expansion for the 
multicenter integral. 

For finite truncation orders ^, Guseinov's rearrangements of the truncations 
( 17.3b are one-range addition theorems since the vectors r and r' are completely sepa- 
rated. Thus, these rearranged truncations can safely be used in multicenter integrals, 
and it is not necessary to split up the integration contour But in the limit ^ -^ °°, 
these rearrangements lose their convenient one-range property. Consequently, their 
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careless use in a multicenter integral without splitting up the integration contour may 
produce either a wrong or a divergent result. 

It makes sense to assume that the ultimate two-range nature of the rearrangements 
of the truncations (17.3b becomes noticeable already for sufficiently large, but finite 
values of jY . It is therefore conceivable that multicenter integrals, whose integrands 
contain a rearrangement of a truncations ( 17.31 ), are semiconvergent with respect to a 
variation of J/ . Obviously, the exact behavior of these integrals as jV -^ °a certainly 
does not only depend on the approximation to the addition theorem, but also on the 
remaining integrand. This certainly makes a detailed analysis even more difficult. 

These considerations are for the moment essentially speculation. A sufficiently 
rigorous analysis of this possible semiconvergence cannot be done yet. We lack some 
necessary mathematical tools such as a two-center analog of the leading order asymp- 
totic expansion ( I7.2l i which was so very useful in the one-center case. 



8 Summary and Outlook 

Infinite dimensional function spaces and in particular Hilbert spaces are of consid- 
erable importance not only in quantum theory, but also in approximation theory and 
in functional analysis. Accordingly, there is an extensive literature discussing these 
function spaces, both from a purely mathematical point of view and also from the 
perspective of quantum physics. 

It is an obvious idea to try to understand the properties of these admittedly com- 
plicated function spaces by emphasizing the analogies to the simpler n-dimensional 
real and complex vector spaces R" and C", respectively. Often, this approach pro- 
vides valuable insight. Nevertheless, there are dangers. Infinite dimensional function 
spaces possess certain peculiarities which do not exist in the case of n-dimensional 
vector spaces. Naive generalizations and overly optimistic analogies can therefore be 
badly misleading. 

For example, in the n-dimensional vector spaces M" or C" any set of n linearly 
independent vectors can be used as a basis, which means that every vector belonging 
to either R" or C" can be expressed as a linear combination of n linearly independent 
vectors. There is also no principal problem if we want to replace one basis by another: 
All linearly independent sets of n vectors are in that respect equivalent. Orthogonal- 
ity of the basis vectors in n-dimensional vector spaces is convenient since it greatly 
simplifies certain operations, but it is not really essential. 

In infinite dimensional function spaces, the situation is much more complicated. 
Firstly, a basis now consists of an infinite number of elements which always raises 
the question of convergence. Moreover, there are different convergence types, which 
are in general incompatible. Secondly, a basis in a function space has to be complete, 
i.e., the span of this basis has to be dense in the function space. 

As discussed in Appendix |C] the completeness of a basis suffices to guarantee 
that finite expansion of the type of ( IC.lb exist. In addition, completeness implies that 
the corresponding mean square deviation ( IC.2b can be made as small as we like by 
increasing the length of the finite expansion. Therefore, it is a seemingly obvious con- 
clusion that the completeness of a basis implies the existence of infinite expansions of 
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the type of ( IC.3l l in terms of this basis. Unfortunately, this conclusion is wrong. The 
existence of an infinite expansion is only guaranteed if the basis is not only complete, 
but also orthogonal. These facts are well known, but nevertheless often ignored. They 
are also the basis of this article. 

In 1I41I - I44 [461, Guseinov derived one-range addition theorems for Slater-type 



functions ;)f^^(j3,r ± r') with integral and nonintegral principal quantum numbers 
A^ by expanding them in terms of his Laguerre-type functions {<;f^"'^(7, J") } ,, ■ For 
a given k — —1,0,1,2,..., these functions defined by (12.13b are complete and or- 
thonormal in the corresponding weighted Hilbert space L\(R^) defined by ( 12.151 ). 
Guseinov's approach is mathematically sound as long as the Slater-type function 
XJ^i^{P,r±r') belongs to L^^. (R^). Such a one-range addition theorem converges in 
the mean with respect to the norm of the Hilbert space L\(R^), but not necessarily 
pointwise. 

However, Guseinov considered it to be advantageous to replace in his A;-dependent 
addition theorems the complete and orthonormal functions {,tf^'"^(7, '")} » , whose 
radial parts are according to ( 12.131 ) essentially generalized Laguerre polynomials, 
by complete, but nonorthogonal Slater-type functions {Zm"^(7i'')} f with integral 
principal quantum numbers via (13.21 1. Ultimately, Guseinov's approach corresponds 
to the replacement of a Laguerre series of the type of (14.4b by a power series of 
the type of (14.5b . Unfortunately, this approach is not necessarily legitimate and can 
easily lead to nontrivial convergence and existence problems, which had already been 
discussed in 1.207.1 . albeit in a less detailed way. 

Because of these principal problems, Guseinov cannot tacitly assume that the 
transformation of his one-range addition theorems, which are expansions in terms 
of his complete and orthonormal Laguerre-type functions, to expansions of the type 
of ( 14.11 ) in terms of the complete, but nonorthogonal Slater-type functions are nec- 
essarily legitimate. This has to be demonstrated explicitly, but Guseinov had only 
done this in a very superficial way. There is considerable evidence that Guseinov and 
co-workers had only done some purely numerical tests. 

The deficiencies of Guseinov's approach become particularly evident in the case 
of the one-center expansion (11.2b for a Slater-type function XnlO^t'') ^i'^^ an in 
general nonintegral principal quantum number A^ in terms of Slater-type functions 
{z^l(7j'")} -r I 1 with integral principal quantum numbers n and an in general dif- 
ferent common scaling parameter 7 7^ j3 > 0. Guseinov |46, Eq. (21)] had constructed 
this expansion by performing the one-center limit of a rearranged truncated addition 
theorem for Slater-type functions |46, Eq. (15)]. This one-center limit was later used 
by Guseinov and Mamedov 1 122] for the construction of series expansions for over- 
lap integrals of Slater-type functions with nonintegral principal quantum numbers. 

It is, however, trivially simple to show that Guseinov's one-center expansion ( 11.21 ) 
does not exist if the principal quantum number A^ is nonintegral. As discussed in Sec- 
tion [T] this follows at once from the simple and yet consequential fact that expan- 
sions in terms of Slater-type functions with integral principal quantum numbers and 
a common scaling parameter are nothing but power series expansions about r = in 
disguise. Moreover, every power series is also a Taylor series for some function (see 
for example 115511 ^ 
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This fact is extremely helpful because the factors, which govern the analyticity 
of a function, are fairly well understood. It is trivially simple to show that the radial 
part of e\p{Yr)XN liP , r) is not analytic in the sense of complex analysis at r = if 
the principal quantum number A^ is not a positive integer satisfying N > L+l. Thus, 
a power series in r for exp{Yr)X]!j ^iP , r) with A^ e M \ N cannot exist, which implies 
that the one-center expansion ( ll.2l i for XNiiP ■ *") 'ioes not exist if A^ G M \ N. 

We arrive at the same conclusion if we apply the transformation formula (14.111 ) 
to the fc-dependent expansions ( 15.4b for XnlH^^'') ^^ terms of Guseinov's complete 
and orthonormal functions or to the equivalent Laguerre series ( 15.11) for z^e"', from 
which (15.4b was derived. As discussed in Section |5] the large index asymptotics of 
the coefficients in the Laguerre series (15.1b clearly shows that the expansion (11.2b for 
XiJii^ ■ *") ™ terms of Slater-type functions does not exist if A^ is nonintegral. 

In principle, the same strategy could also be pursued if one-range addition the- 
orems for Xj^^(j3,r± r'), which are A;-dependent two-center expansions in terms 
of Guseinov's complete and orthonormal function {<:f^,™ (7, '")} » , are rearranged. 
The expansion coefficients of these addition theorems are overlap integrals. On the 
basis of our current level of mathematical understanding, it is, however, very diffi- 
cult or even practically impossible to determine the large index asymptotics of these 
complicated integrals. Fortunately, by means of an indirect approach it is neverthe- 
less possible to arrive at some useful conclusions about the validity of Guseinov's 
rearrangements of his one-range addition theorems. 

My key argument in Sectionals again analyticity in the sense of complex anal- 
ysis. If a function /: C — ;> C has a singularity somewhere in the complex plane, 
its power series about the origin cannot have an infinite radius of convergence. In 
contrast, singularities affect Laguerre expansions only if they are nonintegrable and 
located on the integration contour [0,oo). 

Both exp(— j8|r ±r'|) as well as exp(7r)exp(— j3|r±r'|) are for r' > analytic 
at r = 0. Accordingly, these functions possess power series expansions in r which 
converge in a vicinity of the expansion point r = 0. However, exp(— j3|r ± r'|) is 
singular for r ± r' = 0. Thus, the power series expansions for both exp ( — j3 1 r ± r' | ) 
and exp(7r)exp(— j3|r±r'|) converge only for r < r'. This implies that a one-range 
addition theorem for the Is function exp(— /3|r±r'|), which converges pointwise for 
all r,r' e M? or equivalently for all r,r' e [0,oo), cannot exist. 

These considerations apply also to other Slater-type functions Xnl (i3 , r± /•') with 
integral or nonintegral principal quantum numbers A^. If r' > 0, these functions as well 
as the related functions exp(7r);ifj^^(j3,r±r') are analytic at r = 0, but they also have 
a singularity at r ± r' = 0. Therefore, an expansion of the type of (14.1b in terms of 
Slater-type functions with integral principal quantum numbers and a common scaling 
parameter, that converges pointwise for all r, r' G R^, cannot exist. Instead, we obtain 
the small r part of a two-range addition theorem, either by doing a Taylor expansion 
of Xni il^-J^f') about r == or by applying the transformation formula ( 14.111 ) to the 
angular projections ( 13.101 ). 

The results of Sections |5] and |6] can be summarized as follows: The one-center 
expansion (11.2b for Xnl(P''') 'ioes, not exist if A^ is nonintegral, and the expansion 
(14.1b for ;i;^^(j3,r ± r'), which looks like a one-range addition theorem, has a two- 
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range form since it converges only for \r\ < \r'\. Therefore, one might be tempted to 
dismiss Guseinov's rearrangements of expansions in terms of generalized Laguerre 
polynomials as being both useless and dangerous. 

However, the situation is more complicated than it may look at first sight. As 
discussed in Section |7] Guseinov apparently never rearranged infinite one-center ex- 
pansions of the type of ( 15.41 ) or infinite one-range addition theorems of the type of 
(I3.4l i. although this can be done in a systematic way with the help of the transforma- 
tion formula ( 14.111 ). 

Instead, Guseinov exclusively rearranged^n/fe truncations of his infinite expan- 
sions such as the truncation ( 17.31 1 of the addition theorem (13.7b or the truncation ( 17.11 ) 
of its one-center limit ( I5.4l i. all with a finite truncation order .yV. This is highly con- 
sequential. If we replace in a truncated Laguerre series of the type of ( I4.6l l the gen- 
eralized Laguerre polynomials by powers, for example via ( I4.8l l. there can be no 
convergence or existence problems because the resulting expression is simply a poly- 
nomial. 

Accordingly, Guseinov's rearrangements of truncated one-center and two-center 
expansions produce approximations to XNLiP^'') ^'^'^ XmlH^t''^''')^ respectively, 
that consist of a finite number of Slater-type functions XneiY^ '') with integral princi- 
pal quantum numbers and a common scaling parameter. Since only a finite number 
of terms is transformed in these rearrangements, there can be no existence problems. 
Nevertheless, Guseinov's approach is affected bu some serious problems beyond the 
possible numerical stability problems discussed in Section[3] 

In Section |2] it is shown that rearrangements of the truncations ( 17.1b of the one- 
center expansion ( 15.41 ) are semiconvergent with respect to the truncation order .A' if 
the principal quantum number A^ of the Slater-type function is nonintegral. It follows 
from the leading order asymptotic approximation ( 17.21 ) that the accuracy of these ap- 
proximations first increases with increasing truncation order ,_A^, but for larger values 
of jy the accuracy decreases again, and in the limit of infinite truncation orders, 
everything goes to pieces. 

The situation is not nearly so well understood in the case of rearrangements of 
truncations ( 17.31 ) of the A:-dependent addition theorem ( 13.7b . For finite truncation or- 
ders jy , these rearrangements are one-range addition theorems, but as ,yV — > o° they 
converge to two-range addition theorems. 

It makes sense to assume that the ultimate two-range nature of these rearrange- 
ments becomes noticeable in integrals already for sufficiently large, but finite values 
of ,yf . Thus, it is conceivable that at least certain multicenter integrals containing 
such rearranged truncations might turn out to be semiconvergent with respect to jV . 
For the moment, these considerations are essentially speculations, since substantial 
mathematical knowledge is still lacking. 

As is well known from the literature, the use of semiconvergent expansions can 
offer computational benefits under favorable circumstances. In the case of Guseinov's 
rearranged truncations of one-center and two-center expansions this may also be the 
case. However, the use of semiconvergent expansions or of other approximations, 
whose limit of infinite truncation order either does not exist or has undesirable fea- 
tures, clearly involves some risks. So, before using such an approximation in an ac- 
tual calculation, we firstly must try to understand the inherent risks, and secondly. 
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we must convince ourselves that we will be able to handle these risks. It would be 
extremely negligent to ignore these risks and treat a semiconvergent expansion like a 
convergent expansion. 

The examples mentioned above or the ones given in Appendix |E| should suffice 
to convince even a skeptic reader that the apparent convergence of the leading terms 
of an infinite series alone does not prove anything. We also need some additional 
mathematical insight indicating convergence. 

Numerical demonstrations have obvious limitations. It is always desirable to aug- 
ment them by sufficiently rigorous mathematical investigations. Unfortunately, rig- 
orous proofs are extremely difficult in a research topic as complex and computer 
oriented as electronic structure theory. In the vast majority of all problems in elec- 
tronic structure theory, a rigorous mathematical analysis is out of reach and we have 
to be content with numerical demonstrations, in spite of their obvious limitations and 
their capacity of misleading us. However, the construction of addition theorems and 
their application in multicenter integrals is a mathematical problem, and - as shown 
in this article - quite a few things can be understood on the basis of mathematical 
considerations. Guseinov's strategy of relying entirely on numerical demonstrations 
without trying to understand the underlying mathematics is not acceptable from a 
methodological point of view. 



A General Aspects of Series Expansions 

For the sake of simplicity, let us consider functions F: C — > C A series expansion E^^oMn%(z) for such 
a function F(z) requires a sequence {«„}"^o °^ fixed coefficient and a sequence {'%(z)}^=o °^ known 
functions ^/„: C^ C. 

To make the expansion E"=o""'^i(z) useful, it has to represent F{z) in some sense. Thus, we assume 
that Y.n^oUii'^n(z) converges to F(z) according to some specified convergence type and we write: 

F(z) = Y.u„^„{z). (A.l) 

11=0 

Numerous different convergence types occur in practice. Important examples are pointwise convergence, 
which is typical of classical complex analysis or also of two-range addition theorems, convergence in 
the mean with respect to the norm of some Hilbert space, which is typical of most one-range addition 
theorems, or even distributional or weak convergence in the sense of Schwartz 1 174]. 

Therefore, the indiscriminate use of the "=" sign for all convergence types is potentially misleading 
since it suggests a uniqueness which does not exist. Depending on the convergence type, an "=" sign 
can have a completely different meaning, or to put it differently, series expansions of the type of jA. It 
can have very different mathematical properties. Moreover, different convergence types are in general 
incompatible, i.e., the convergence of a series expansion for a given convergence type does not imply that 
this series converges also with respect to another convergence type. 

It is, however, overly restrictive to insist that a series expansion of the type of jA. It must converge in 
some sense to be practically useful. Divergent, but summable series expansions are simply to useful to be 
discarded, not only in mathematics, but in particular in quantum physics (see for example 1 150] and ref- 
erences therein). In numerous scientific applications, there is no alternative to the summation of divergent 
series. A condensed review of divergent series and their summation can be found in |207, Appendices A 
andB]. 

Let us now assume that we have a sequence of approximate expressions of the following kind: 



JV 

I 

n=0 



Fn (z) = £ "f %Az), W e No . (A.2) 
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The superscript N in u„ indicates that the coefficients may depend explicitly on the summation limit N, 

i.e., we in general have u„ y^ ul, =^ Un =^ ■■■ for fixed n,N S No. 

Let us now also assume that the approximants Fjv(z) converge to F(z) as W — > t» in some sense. This 
raises the question whether the resulting expression F[z) = limAi^oo iw (z) constitutes an expansion of F(z) 
in terms of the functions {"2^1(2) },7=o- 

The answer is that this is in general not true. The convergence of the approximants /vv(z) to F(z) 
as W — > 00 only means that we can make the difference between Fn(z) and F(z), whose exact meaning 
depends on the convergence type, as small as we like by increasing N as much as necessary. This does not 
guarantee that the coefficients u„ in )A.2) have for all « e No unique limits u„ = ii„ = limAi^oo «,, . 
Thus, the existence of a convergent sequence of approximants of the type of )A.2t does not imply the 
existence of a series expansion of the type of jA.lt . 



B Glory and Misery of Power Series 

It is probably justified to claim that power series 

/W = E ^-^(^■-zoV = t '^Az-zoY (B.l) 

v=0 ^ ■ v=0 

are the most important analytical tools not only in mathematical analysis, but also in the mathematical 
treatment of scientific and engineering problems. 

As is well known, power series can be differentiated term-by-term under relatively mild conditions. 
Accordingly, it is an obvious idea to try to solve differential equations in terms of power series. A large 
part of special function theory consists of the construction and analysis of power series solutions to the 
ordinary differential equations that are of relevance in mathematical physics. Since power series can also 
be integrated term-by-tenn under very mild conditions, they are also indispensable for the construction of 
explicit expressions for integrals involving functions that possess power series representations. 

While the usefulness of power series in analytical manipulations cannot be overemphasized, it is 
nevertheless also true that a power series representation for a given function is at best a mixed blessing 
from a purely numerical point of view. The problem is that power series expansions converge in circles 
about the expansion point. As is well known, the radius of such a circle of convergence is determined by 
the location the closest singularity of the function under consideration. Thus, the radius of convergence of 
a power series expansion can be zero, finite and infinite. 

In general, a power series in z~za is numerically useful only if z and zo differ slightly, i.e., in the 
immediate vicinity of the expansion point zo- Then, a few terms of the series normally suffice to produce 
excellent approximations. But close to the boundary of its circle of convergence, the rate of convergence 
of such a power series expansion can become prohibitively slow. In my opinion, the current popularity of 
Pade approximants, which normally converge much more rapidly than the partials sums from which they 
are constructed (see for example the monograph of Baker and Graves-Morris |5] and references therein), is 
largely due to the combined effect of the undeniable analytical usefulness of power series and their (very) 
limited usefulness as computational tools. 

A nonzero, but finite radius of convergence of a power series can also cause serious problems in 
integrals. As is well known, a series expansion for the integrand can be integrated termwise if it converges 
uniformly for the whole integration interval. Otherwise, we have to be prepared that termwise integration 
either generates a wrong result or even a divergent series expansion for the integral. Thus, if we want to 
replace a part of an integrand by its power series expansion, we are on the safe side only if the integration 
interval is completely contained in the circle of convergence. In particular in the case of infinite or semi- 
infinite integration intervals, this is normally not the case. 



C Orthogonal Expansions 

The analyticity of a function / : C — > C in the sense of complex analysis is undeniably a highly desirable 
feature. In the interior of suitable subsets of C, / can be represented by power series expansions which 
converge pointwise and uniformly, and it is also comparatively easy to compute derivatives of / in this way. 
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Neveitheless, it is often advantageous to use instead of power series alternative expansions that converge 
in a weaJcer sense. 

Let "y be an infinite dimensional vector space with inner product (■!■): T^ x y — > C and its corre- 
sponding norm || ■ || : "Z' — > 1R+ defined by ||/|| = (/|/)''^ with / e "V . If every Cauchy sequence in Y 
converges with respect to the norm || ■ || to an element of "/', then "V is called a Hilbert space. 

Let us now assume that / is an element of some Hilbert space ,J^ and that the functions { (p,„ D^^q are 
linearly independent and complete in Jif. Then, we can construct approximations 

M 

„(M) 



Y.o"'(p„ (C.l) 



to /, where M is a finite integer. The expansion coefficients C,„ , which in general depend on the summa- 
tion limit M, are chosen in such a way that the mean square deviation 

II/-/m||- = if-fM\f-fM) (C.2) 

becomes minimal. 

The finite approximation )C. It converges to / as M — > oo if )C.2t can be made as small as we like by 
increasing M. It therefore looks natural to assume that / possesses an infinite expansion 

/ = £ C„,(p,„ (C.3) 

in terms of the linearly independent and complete functions {<Pm}'^,^Q with well defined expansion coeffi- 
cients C„, = HmM^oo Cm . 

As discussed in Appendix|A] this naturally looking assumption is not necessarily true. In general, the 
coefficients C,„ in )C. U do not only depend on m, /, and {<Pm}Z=0' '"^' ^'*'° "" ^^^ summation limit M. It 
is not a priori clear whether the coefficients C„, in jC.U possess well defined limits C,„ = limM^oo C„, , 
or whether an infinite expansion of the type of )C.3t exists. In fact, expansions of that kind may or may 
not exist. 

It is one of the central results of approximation theory that for arbitrary / 6 J^ the mean square 
deviation ||/ — /m|P becomes minimal if the functions {<Pm}Z=o ^^ ^'^^ °^y lii^^arly independent and 
complete, but also orthonormal satisfying {(p,„\(p„i) = S^„i for all m,m' 6 Mq, and if the coefficients are 
chosen according to C,„ = {%,]/) (see for example |24, Theorem 9 on p. 51]). 

If the functions {<Pm}'^,^o are complete and orthonormal in the Hilbert space ,#' and if the expansion 

coefficients in )C. it are chosen according to Cm — (<Pm|/), then the expansion coefficients do not depend 
on M. Thus, the limit M — > <» is possible, and / G J^ possesses an infinite series expansion 

/= i;(<P».|/)'Pm (C.4) 

m=0 

and this expansion converges in the mean with respect to the norm || ■ || of Jif. 

The fact, that the completeness of a function set {fPmJm^o ™ ^^ infinite dimensional Hilbert space ,^ 
alone does not suffice to guarantee the existence of expansions of the type of jC.3t . is highly consequential. 
Nevertheless, it is occasionally overlooked, although this insufficiency is well documented both in the 
mathematical literature (see for example |24, Theorem 10 on p. 54] or 1137, Section 1.4]) as well as 
in the literature on electronic structure calculations 1 143- 148]). Horrifying examples of nonorthogonal 
expansions with pathological properties can be found in (144, Section III. I]. 

If the Hilbert space j*f is an infinite dimensional vector space consisting of function f,g: C — > C, 
the inner product {f\g) of ,#" is usually identified with an integral J^ w(z)[f{!)]* g{z)dz, where w{z) is an 
appropriate positive weight function. Moreover, the complete orthonormal functions {(Pm}Z=o '^ '^ ^^ 
normally related to a suitable subclass {.i'^m(z)})^=o °t'the classical orthogonal polynomials of mathemat- 
ical physics, as specified by the integration limits a and b and the weight function w(z). In this case, the 
general orthogonal expansion )C.4t boils down to the expansion of a function /(z) in terms of orthogonal 
polynomials: 

/(Z) = I 'gm^^Az) (C.5) 

/ii=0 
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It is generally accepted that orthogonal expansions are extremely useful mathematical tools and that 
they have many highly advantageous features. This is, however, not the whole truth, in particular if we want 
to approximate functions. Hilbert space theory only guarantees that an orthogonal expansion converges in 
the mean with respect to the corresponding norm || ■ ||, but not necessarily pointwise or even uniformly. 
Thus, convergence in the mean is a comparatively weak form of convergence, and orthogonal expansions 
are not necessarily a good choice if we are predominantly interested in the local properties of a function. 
However, convergence in the mean is usually completely satisfactory for the evaluation of integrals. 



D Generalized Laguerre Polynomials 



The surface spherical harmonics Y"'{6,(j)) are complete and orthonormal with respect to an in tegration 
over the surface of the unit sphere in M? (an explicit proof can for instance be found in \lm„ Section 
111.7.6]). Since more complex Hilbert spaces can be constructed by forming tensor products of simpler 
Hilbert spaces (see for example f 164, Section 11.6.5]), we only have to find suitable radial functions that 
are complete and orthogonal with respect to an integration from to oo (see also 1147, Lemma 6 on p. 31]). 
Thus, we more or less automatically arrive at function sets based on the generalized Laguerre polynomials. 
The generalized Laguerre polynomials L„ (z) with 9?(a) > — 1 and n e No are orthogonal with 
respect to an integration over the positive real semiaxis [0,°°) with weight function w{z) = z" exp(— z). In 
the mathematical literature, they are defined either via the Rodrigues relationship 



n\ az 



;,"+«] 



(D.l) 



or as a terminating confluent hypergeometric series i Fi 

l\ 



(«).,) ^ (f^,f,(_„;« + l;,) = (f^ f ±^l_ 



(D.2) 



"' v=0 («+l)v v! 
The generalized Laguerre polynomials satisfy for 9?(a) > — 1 and m, n G No the orthogonality relationship 



i: 



'.«e-Lif'(z)Li«'(z)d.= £(^±f±ll5„ 



(D.3) 



Accordingly, the polynomials 



^!^\z) = 



r{a + n+\] 



1/2 



4"'W, 



n e No , a > - 1 , 



(D.4) 



are for 9t(a) > — 1 orthonormal with respect to an integration over the interval [0,oo) involving the weight 
function z" exp( — z) : 



t-^z''J^L"\z)J^li"\z)dz = S,„ 



(D.5) 



Alternatively, we can also use the functions 



{c)t^^ _ 



^^iz) 



1/2 



e-^/2.«/2Lr)(z), neNo, 



_r(a + «+l)_ 
which are for 9t(a) > — 1 orthonormal with respect to an integration over the interval [0, oo): 



I 



<pi"\z)1'i"\z)dz = S„, 



(D.6) 



(D.7) 



The completeness of the generalized Laguerre polynomials in the weighted Hilbert space 



lU^ 



([0,00)) = {/: [0,oo)^c\J\-'z"\f{z)\^dz<oo] 



(D.i 



is a classic result of mathematical analysis (see for example Il37l p. 33], Il72l pp. 349 - 351], or [ 1 89 , pp. 
235 - 238]). 
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In general, LagueiTe expansions converge only in the mean, but not necessarily pointwise (see for 
example [^). Additional conditions, which a function has t o sati sfy in order to guarantee that its Laguerre 
expansion converges pointwise, were discussed by Szego llSd Theorem 9.1.5 on p. 246] (see also I29l 
Appendix]). 

In this article, exclusively the mathematical notation is used. A different convention for Laguerre 
poly nomials is frequently used in the quantum mechanical literature. For example, Bethe and Salpeter 
rni Eq. (3.5)] define associated Laguerre functions [^"/(zJJrc with n,m G No via the Rodrigues-type 
relationships 



d" 

d?" 



[CW]bs = x;;r[^«W]Bs. (D-9a) 



[Uz)]ss=^'^K'z"\- (D.9b) 



dz« 



Comparison of jD.U and jD.9b) implies: 



The convention of Bethe and Salpeter 1 11] is also used in the books by Condon and Shortley ^211, Eqs. (6) 
and (9) on p. 115] and by Condon and Odaba§i |20, Eq. (2) on p. 189] as well as in the numerous articles 
by Guseinov and his coworkers. 

In my opinion, the use of associated Laguen'e functions defined by jD.9) is not recommendable. It 
follows from JD.IQ) that these functions cannot express generalized Laguen'e polynomials L„ with non- 
integral superscripts a. This is both artificial and unnecessary. For example, the eigenfunctions Q'^'^{ji,r) 
of the Hamiltonian /3""V- — ji^r^ of the three-dimensional isotropic harmonic osci llator contain gener- 
alized Laguerre polynomials in r^ with half-integral superscripts (see for example [194, Eq. (5.4)] and 
references therein). 



E Semiconvergence 

If the sum of the leading terms of an infinite series seem to approach the connect limit, it looks natural to 
conclude that increasing the number of terms will improve the accuracy of the approximation, and that 
ultimately the partial sums will converge to the correct limit. Unfortunately, this assumption is overly 
optimistic and in spite of its apparent plausibility not necessarily fitie. 

Many series are known whose partial sums initially seem to converge. If however, further terms are 
included, the accuracy decreases, and ultimately the sequence of partial sums diverges. In th e liter ature, 
this phenomenon is well established and usually called semiconvergence (see for exa mple 11621 p. 2, 
Footnote^]). To the best of my knowledge, this terminology was introduced by Stieltjes Il85ll already in 
1886. 

Semiconvergence is best known in connection with factorially divergent asymptotic inverse power 
series for special functions. In Arfken's book (j Chapter 5.10], one can find a comprehensive disc ussio n of 
the semiconvergence of the divergent asymptotic series of the incomplete gamma functi on F (a,z) ll63l Eq. 
(8.11.2)] and its special case, the asymptotic series for the exponential integral £1 (z) fl53', Eq. (6.12.1)]. 
In my opinion, Arfken's treatment is well suited as a first introduction to this topic. Other examples of 
semiconvergent series are the factorially divergent asymptotic series for the complementary error function 
erfc(z) (163, Eq. (7.12.1)], and the modified Bessel and Whittaker functions of the second kind Kv{z) and 
Wk,m(z). respectively 1 163, Eqs. (10.40.4) and (13.19.3)]. 

If the argument of such a divergent asymptotic inverse power series is sufficiently large, then the trun- 
cation of such a divergent series in the vicinity of the minimal term can lead to excellent approximations 
to the function it represents. Nevertheless, these partial sums diverge if further terms beyond the minimal 
term are included (see for example |34, Figure 2.2 on p. 35]). 

Accordingly, the accuracy, which can be obtained by truncating a semiconvergent inverse power series 
in the vicinity of the minimal term, depends crucially on the magnitude of the argument. If the argument is 
large, excellent approximations can often be obtained. The situation is not so good if the argument is small, 
because then the truncation of a semiconvergent series can only provide relatively crude approximations. 
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But even for small ai'guments, it is often possible to obtain very accurate approximations by using the 
partial sums of a semiconvergent series as input data in nonlinear sequence transformations fl98l,l212l]. 

Semiconvergent series occur also quite abundantly in quantum mechanical perturbation expansions. 
For example, Ahlrichs [2j] showed that the total energy of interacting molecular systems A and B can be 
expressed by a semiconvergent series, the so called l//?-expansion. 

Other phenomena closely resembling semiconvergence are also known. Baumel, Crocker, and Nuttall 
131 showed that in the case of scattering calculations with complex basis functions low order approxima- 
tions can produce good results although the whole scheme ultimately diverges, and Gautschi 1 32] observed 
initial apparent convergence to the wrong limit in the case of continued fractions for Kummer functions. 



References 

1. Abramowitz, M., Stegun, I. A. (eds.): Handbook of Mathematical Functions. 
National Bureau of Standards, Washington, D. C. (1972) 

2. Ahlrichs, R.: Convergence properties of the intermolecular force series {l/R- 
expansion). Theor. Chim. Acta 41, 7-15 (1976) 

3. Arfken, G.B.: Mathematical Methods for Physicists, 3 edn. Academic Press, 
Orlando (1985) 

4. Askey, R., Wainger, S.: Mean convergecne of expansions in Laguerre and Her- 
mite series. Amer J. Math. 87, 695 - 708 (1965) 

5. Baker Jr, G.A., Graves-Morris, P.: Pade Approximants, 2 edn. Cambridge U. 
P, Cambridge (1996) 

6. Barnett, M.P, Coulson, C. A. : The evaluation of integrals occurring in the theory 
of molecular structure. Part I & II. Phil. Trans. Roy. Soc. London A 243, 221 - 
249(1951) 

7. Baumel, R.T., Crocker, M.C., Nuttall, J.: Limitations of the method of complex 
basis functions. Phys. Rev. A 12, 486 - 492 (1975) 

8. Beckermann, B., Kalyagin, V., Matos, A.C., Wielonsky, E: How well does the 
Hermite-Pade approximation smooth the Gibbs phenomenon? Math. Comput. 
80,931-958(2011) 

9. Beckermann, B., Matos, A.C., Wielonsky, R: Reduction of the Gibbs phe- 
nomenon for smooth functions with jumps by the e-algorithm. J. Comput. Appl. 
Math. 219, 329 - 349 (2008) 

10. Bender, CM., Weniger, E.J.: Numerical evidence that the perturbation expan- 
sion for a non-Hermitian ^i^-symmetric Hamiltonian is Stieltjes. J. Math. 
Phys. 42,2167-2183(2001) 

11. Bethe, H.A., Salpeter, E.A.: Quantum Mechanics of One- and Two-Electron 
Atoms. Plenum Press, New York (1977). Originally published by Springer- 
Verlag in Handbuch der Physik Vol. XXXV, Atome I, BerHn (1957). 

12. Biedenharn, L.C., Louck, J.D.: Angular Momentum in Quantum Physics. 
Addison-Wesley, Reading, Mass. (1981) 

13. Borghi, R.: Evaluation of diffraction catastrophes by using Weniger transfor- 
mation. Opt. Lett. 32, 226 - 228 (2007) 

14. Borghi, R.: Joint use of the Weniger transformation and hyperasymptotics for 
accurate asymptotic evaluations of a class of saddle-point integrals. Phys. Rev. 
E 78, 026,703-1 - 026,703-1 1 (2008) 



54 



15. Borghi, R.: Joint use of the Weniger transformation and hyperasymptotics for 
accurate asymptotic evaluations of a class of saddle-point integrals. 11. Higher- 
order transformations. Phys. Rev. E 80, 016,704-1 - 016,704-15 (2009) 

16. Brezinski, C: Extrapolation algorithms for filtering series of functions, and 
treating the Gibbs phenomenon. Numer. Algor 36, 309 - 329 (2004) 

17. Brezinski, C, Redivo Zaglia, M.: A review of vector convergence acceleration 
methods, with applications to linear algebra problems. Int. J. Quantum Chem. 
109,1631-1639(2009) 

18. Browne, J.C.: Molecular wave functions: Calculation and use in atomic and 
molecular processes. Adv. At. Mol. Phys. 7, 47 - 95 (1971) 

19. Cizek, J., Zamastil, J., Skala, L.: New summation technique for rapidly diver- 
gent perturbation series. Hydrogen atom in magnetic field. J. Math. Phys. 44, 
962 - 968 (2003) 

20. Condon, E.U., Odabagi, H.: Atomic Structure. Cambridge U. P., Cambridge 
(1980) 

21. Condon, E.U., Shortley, G.H.: The Theory of Atomic Spectra. Cambridge U. 
P, Cambridge (1970) 

22. Cvetic, G., Yu, J.Y.: Borel-Pade vs Borel- Weniger method: A QED and a QCD 
example. Mod. Phys. Lett. A 15, 1227 - 1235 (2000) 

23. Dalgarno, A.: Integrals occurring in problems of molecular structure. Math. 
Tables Aids Comput. 8, 203 - 212 (1954) 

24. Davis, H.F.: Fourier Series and Orthogonal Functions. Dover, New York (1989). 
Originally published by Allyn and Bacon, Boston (1963) 

25. Doha, E.H.: On the connection coefficients and recurrence relations arising 
from expansions in series of Laguerre polynomials. J. Phys. A 36, 5449 - 5462 
(2003) 

26. Erdelyi, A., Magnus, W., Oberhettinger, F, Tricomi, E.G.: Higher Transcenden- 
tal Functions, vol. 11. McGraw-Hill, New York (1953) 

27. Filter, E.: Analytische Methoden zur Auswertung von Mehrzentren- 
Matrixelementen in der Theorie der Molekiilorbitale bei Verwendung expo- 
nentialartiger Basissatze. Ph.D. thesis, Fachbereich Chemie und Pharmazie, 
Universitat Regensburg (1978) 

28. Filter, E., Steinborn, E.O.: Extremely compact formulas for molecular one- 
electron integrals and Coulomb integrals over Slater-type atomic orbitals. Phys. 
Rev. A 18, 1-11 (1978) 

29. Filter, E., Steinborn, E.O.: A matrix representation of the translation operator 
with respect to a basis of exponentially declining functions. J. Math. Phys. 21, 
2725-2736(1980) 

30. Fock, v.: Naherungsmethoden zur Losung des quantenmechanischen 
Mehrkorperproblems. Z. Physik 61, 126 - 148 (1930) 

31. Gaunt, J.A.: The triplets of helium. Phil. Trans. Roy. Soc. A 228, 151 - 196 
(1929) 

32. Gautschi, W. : Anomalous convergence of a continued fraction for ratios of kum- 
mer functions. Math. Comput. 31, 994 - 999 (1977) 

33. Gautschi, W.: Is the recurrence relation for orthogonal polynomials always sta- 
ble? BIT 33, 277 - 284 (1993) 



55 

34. Gil, A., Segura, J., Temme, N.M.: Numerical Methods for Special Functions. 
SIAM, Philadelphia (2007) 

35. Gottlieb, D., Orszag, S.A.: Numerical Analysis of Spectral Methods: Theory 
and Applications. Regional Conference Series in Applied Mathematics. SIAM, 
Philadelphia (1977) 

36. Gradshteyn, I.S., Rhyzhik, I.M.: Table of Integrals, Series, and Products, 5th 
edition edn. Academic Press, Boston (1994) 

37. Grecchi, V., Maioli, M., Martinez, A.: Pade summability of the cubic oscillator. 
J. Phys. A 42, 425,208-1 - 425,208-17 (2009) 

38. Grotendorst, J., Weniger, E.J., Steinborn, E.O.: Efficient evaluation of infinite- 
series representations for overlap, two-center nuclear attraction, and Coulomb 
integrals using nonlinear convergence accelerators. Phys. Rev. A 33, 3706 - 
3726 (1986) 

39. Guilpin, C, Gacougnolle, J., Simon, Y.: The e-algorithm allows to detect Dirac 
delta functions. Appl. Numer. Math. 48, 27 - 40 (2004) 

40. Guseinov, I., Aydm, R., Mamedov, B.: Computation of multicenter overlap in- 
tegrals with Slater-type orbitals using ^^("'-ETOs. J. Mol. Model. 9, 325 - 328 
(2003) 

41. Guseinov, I.I.: Analytical evaluation of three- and four-center electron- 
repulsion integrals for Slater-type orbitals. J. Chem. Phys. 69, 4990 - 4994 
(1978) 

42. Guseinov, I.I.: Expansion of Slater-type orbitals about a new origin and analyt- 
ical evaluation of multicenter electron-repulsion integrals. Phy. Rev. A 22, 369 
-371(1980) 

43. Guseinov, I.I.: Expansion of Slater-type orbitals about a displaced center and the 
evaluationof multicenter electron-repulsion integrals. Phys. Rev. A 31, 2851 - 
2853(1985) 

44. Guseinov, I.I. : Evaluation of expansion coefficients for translation of Slater-type 
orbitals using complete orthonormal sets of exponential-type functions. Int. J. 
Quantum Chem. 81, 126 - 129 (2001) 

45. Guseinov, I.I.: Computation of molecular integrals over Slater-type orbitals. 
IX. Calculation of multicenter multielectron molecular integrals with integer 
and noninteger n Slater orbitals using complete orthonormal sets of exponential 
functions. J. Mol. Struc. (Theochem) 593, 65 - 69 (2002) 

46. Guseinov, I.I.: New complete orthonormal sets of exponential-type orbitals and 
their application to translation of Slater orbitals. Int. J. Quantum Chem. 90, 1 14 
-118(2002) 

47. Guseinov, 1. 1.: Unified analytical treatment of one-electron multicenter integrals 
of central and noncentral potentials over Slater orbitals. Int. J. Quantum Chem. 
90, 980 - 985 (2002) 

48. Guseinov, I.I.: Addition and expansion theorems for complete orthonormal sets 
of exponential-type orbitals in coordinate and momentum representations. J. 
Mol. Model. 9, 135 - 141 (2003) 

49. Guseinov, I.I.: Addition theorems for Slater-type orbitals and their application 
to multicenter multielectron integrals of central and noncentral interaction po- 
tentials. J. Mol. Model. 9, 190 - 194 (2003) 



56 

50. Guseinov, I.I.: Comment on "Evaluation of two-center overlap and nuclear- 
attraction integrals over Slater-type orbitals with integer and noninteger prin- 
cipal quantum numbers". Int. J. Quantum Chem. 91, 62 - 64 (2003) 

51. Guseinov, I.I.: Unified analytical treatment of multicenter multielectron inte- 
grals of central and noncentral interaction potentials over Slater orbitals using 
f "-ETOs. J. Chem. Phys. 119, 4614-4619 (2003) 

52. Guseinov, I.I.: Unified treatment of integer and noninteger n multicenter multi- 
electron molecular integrals using complete orthonormal sets of f "-ETOs. J. 
Mol. Struc. (Theochem) 625, 221 - 225 (2003) 

53. Guseinov, I.I.: Analytical evaluation of multicenter multielectron integrals of 
central and noncentral interaction potentials over Slater orbitals using overlap 
integrals and auxiliary functions. J. Math. Chem. 36, 83-91 (2004) 

54. Guseinov, 1. 1.: Comment on "Calculation of two-center nuclear attraction inte- 
grals over integer and noninteger «-Slater-type orbitals in nonlined-up coordi- 
nate systems". J. Math. Chem. 36, 123 - 127 (2004) 

55. Guseinov, I.I.: Comment on "Evaluation of multicenter electric multipole mo- 
ment integrals over integer and noninteger n-STOs". J. Chin. Chem. Soc. 51, 
1077-1078(2004) 

56. Guseinov, I.I.: Erratum: Use of f "-ETOs in the unified treatment of elec- 
tronic attraction, electric field and electric field gradient multicenter integrals 
of screened Coulomb potentials over Slater orbitals. J. Mol. Model. 10, 233 
(2004) 

57. Guseinov, I.I.: One-range addition theorems for derivatives of Slater-type or- 
bitals. J. Mol. Model. 10, 212 - 215 (2004) 

58. Guseinov, I.I.: Unified analytical treatment of multicentre electron attraction, 
electric field and electric field gradient integrals over Slater orbitals. J. Phys. A 
37, 957 - 964 (2004) 

59. Guseinov, I.I.: Unified analytical treatment of two-electron multicenter inte- 
grals of central and noncentral interaction potentials over Slater orbitals. Int. J. 
Quantum Chem. 100, 206 - 207 (2004) 

60. Guseinov, 1. 1.: Unified treatment of electronic attraction, electric field, and 
electric-field gradient multicenter integrals of screened and nonscreened 
Coulomb potentials using overlap integrals for Slater orbitals. Can. J. Phys. 
82,819-825(2004) 

6 1 . Guseinov, 1. 1.: Unified treatment of multicenter integrals of integer and noninte- 
ger u Yukawa-type screened Coulomb type potentials and their derivatives over 
Slater orbitals. J. Chem. Phys. 120, 9454 - 9457 (2004) 

62. Guseinov, 1. 1.: Use of f "-ETOs in the unified treatment of electronic attrac- 
tion, electric field and electric field gradient multicenter integrals of screened 
Coulomb potentials over Slater orbitals. J. Mol. Model. 10, 19 - 24 (2004) 

63. Guseinov, I.I.: Addition theorems for Slater-type orbitals in momentum space 
and their application to three-center overlap integrals. J. Mol. Model. 11, 124 - 
127 (2005) 

64. Guseinov, 1. 1.: Corrigendum to "One-range addition theorems for derivatives of 
complete orthonormal sets of "F^-ETOs" [J. Mol. Struct.: THEOCHEM 719 
(2005) 53 - 55]. J. Mol. Struc. (Theochem) 726, 285 (2005) 



57 

65. Guseinov, 1. 1.: Evaluation of two- and three-center overlap integrals over com- 
plete orthonormal sets of f "-ETOs using their addition theorems. J. Math. 
Chem. 38,489-493(2005) 

66. Guseinov, I.I.: One-range addition theorems for Coulomb interaction potential 
and its derivatives. Chem. Phys. 309, 209 - 213 (2005) 

67. Guseinov, 1. 1.: One-range addition theorems for derivatives of complete or- 
thonormal sets of f "-ETOs. J. Mol. Struc. (Theochem) 719, 53 - 55 (2005) 

68. Guseinov, I.I.: One-range addition theorems for derivatives of integer and non- 
integer M Coulomb Yukawa type central and noncentral potentials and their 
application to multicenter integrals of integer and noninteger n Slater orbitals. 
J. Mol. Struc. (Theochem) 757, 165 - 169 (2005) 

69. Guseinov, 1. 1.: One-range addition theorems for Yukawa-like central and non- 
central interaction potentials and their derivatives. Bull. Chem. Soc. Japan 78, 
611-614(2005) 

70. Guseinov, 1. 1.: New complete orthonormal sets of hyperspherical harmonics and 
their one-range addition and expansion theorems. J. Mol. Model. 12, 757 - 761 
(2006) 

71. Guseinov, I.I.: One-range addition theorems for combined Coulomb and 
Yukawa like central and noncentral interaction potentials and their derivatives. 
J. Math. Chem. 39, 253 - 258 (2006) 

72. Guseinov, 1. 1.: Expansion formulae for two-center integer and noninteger n STO 
charge densities and their use in evaluation of multi-center integrals. J. Math. 
Chem. 42,415-422(2007) 

73. Guseinov, 1. 1.: Unified treatment of complete orthonormal sets for wave func- 
tions, and Slater orbitals of particles with arbitrary spin in coordinate, momen- 
tum and four-dimensional spaces. Phys. Lett. A 372, 44 - 48 (2007) 

74. Guseinov, 1. 1.: Unified treatment of complete orthonormal sets of functions in 
coordinate, momentum and four-dimensional spaces and their expansion and 
one-range addition theorems. J. Math. Chem. 42, 991 - 1001 (2007) 

75. Guseinov, I.I.: Combined theory of two-electron nonrelativistic and quasirela- 
tivistic multicenter integrals over integer and noninteger n Slater type orbitals 
using auxiliary functions Qns and Q-n.s- Commun. Math. Comput. Chem. 
(MATCH) 52, 385 - 394 (2008) 

76. Guseinov, I.I.: On the unified treatment of complete orthonormal sets of func- 
tions in coordinate, momentum, and four dimensional spaces and their expan- 
sion and one-range addition theorems. J. Math. Chem. 43, 1024 - 1031 (2008) 

77. Guseinov, I.I.: One-range addition theorems in terms of f "-ETOs for STOs 
and Coulomb- Yukawa like correlated interaction potentials of integer and non- 
integer indices. Chin. Phys. Lett. 25, 4240 - 4243 (2008) 

78. Guseinov, I.I.: Unified treatment of complete orthonormal sets for exponential 
type vector orbitals of particles with spin 1 in coordinate, momentum and four- 
dimensional spaces. J. Math. Chem. 44, 197 - 205 (2008) 

79. Guseinov, I.I.: Unified treatment of complete orthonormal sets of nonrelativis- 
tic, quasirelativistic and relativistic sets of spinor wave functions, and Slater 
spinor orbitals in coordinate, momentum and four-dimensional spaces. J. Math. 
Chem. 44, 839 - 848 (2008) 



58 

80. Guseinov, 1. 1.: Unsymmetrical and symmetrical one-range addition theorems 
for Slater type orbitals and Coulomb-Yukawa-like correlated interaction poten- 
tials of integer and noninteger indices. J. Theor Comput. Chem. 7, 257 - 262 
(2008) 

81. Guseinov, I.I.: Use of Cartesian coordinates in evaluation of multicenter multi- 
electron integrals over Slater type orbitals and their derivatives. J. Math. Chem. 
43, 427 - 434 (2008) 

82. Guseinov, I.I.: Combined theory of nonrelativistic and quasirelativistic atomic 
integrals over integer and noninteger n-Slater-type orbital. J. Theor Comput. 
Chem. 8, 47 - 56 (2009) 

83. Guseinov, I.I.: Errata: "Unsymmetrical and symmetrical one-range addition the- 
orems for Slater type orbitals and Coulomb-Yukawa-like correlated interaction 
potentials of integer and noninteger indices" [Journal of Theoretical and Com- 
putational Chemistry, Vol. 7, No. 2 (2008) 257 - 262]. J. Theor Comput. Chem. 
8, 183 (2009) 

84. Guseinov, I.I.: Evaluation of multicenter multielectron integrals using one- 
range addition theorems in terms of STOs for STOs and Coulomb-Yukawa like 
correlated interaction potentials with integer and noninteger indices. Bull. Ko- 
rean Chem. Soc. 30, 1539 - 1542 (2009) 

85. Guseinov, I.I.: Expansion formulae for one- and two-center charge densities 
over complete orthonormal sets of exponential type orbitals and their use in 
evaluation of multicenter-multielectron integrals. J. Theor. Comput. Chem. 8, 
597 - 602 (2009) 

86. Guseinov, I.I.: Theory of complete orthonormal relativistic vector wave func- 
tion sets and Slater type relativistic vector orbitals in coordinate, momentum 
and four-dimensional spaces. J. Math. Chem. 45, 946 - 952 (2009) 

87. Guseinov, I.I.: Theory of complete orthonormal sets of relativistic tensor wave 
functions and Slater tensor orbitals of particles with arbitrary spin in position, 
momentum and four-dimensional spaces. Phys. Lett. A373, 2178-2181 (2009) 

88. Guseinov, I.I.: Theory of complete orthonormal sets of relativistic tensor wave 
functions and Slater tensor orbitals of particles with arbitrary spin in position, 
momentum and four-dimensional spaces. Phys. Lett. A 373, 2178-2181 (2009) 

89. Guseinov, I.I.: Unified treatment of nonrelativistic and quasirelativistic atomic 
integrals over complete orthonormal sets of f "-exponential type orbitals. J. 
Math. Chem. 45, 1032 - 1040 (2009) 

90. Guseinov, I.I.: Use of auxiliary functions Qlis and Gl^, in evaluation of mul- 
ticenter integrals over integer and noninteger n-Slater type orbitals arising in 
Hartree-Fock-Roothaan equations for molecules. J. Math. Chem. 45, 974 - 980 
(2009) 

91. Guseinov, 1. 1.: Use of Coulomb- Yukawa like correlated interaction potentials of 
integer and noninteger indices and one-range addition theorems for f "-ETO in 
evaluation of potential of electric field produced by molecule. Bull. Korean 
Chem. Soc. 30, 2617 - 2620 (2009) 

92. Guseinov, I.I.: Combined theory of two- and four-component complete or- 
thonormal sets of spinor wave functions and Slater type spinor orbitals in posi- 
tion, momentum and four-dimensional spaces. J. Math. Chem. 47, 391 - 402 



59 

(2010) 

93. Guseinov, 1. 1.: Evaluation of intermolecular interaction energy using one-range 
addition theorems for f^-ETO and Coulomb-Yukawa like correlated interac- 
tion potentials with integer and noninteger indices. J. Math. Chem. 47, 1240 - 
1247 (2010) 

94. Guseinov, 1. 1.: Expansion formulae for two-center charge densities of integer 
and noninteger n generalized exponential type orbitals applied to evaluation of 
multicenter multielectron integrals. J. Math. Chem. 47, 384 - 390 (2010) 

95. Guseinov, 1.1.: Expansion formulae for two-center charge densities of integer 
and noninteger n generalized exponential type orbitals with hyperbolic cosine 
and their use in evaluation of multicenter multielectron integrals. J. Math. 
Chem. 47, 1240-1247 (2010) 

96. Guseinov, 1.1.: Evaluation of potential of electric field produced by molecule 
using symmetrical one-range addition theorems for Coulomb-Yukawa like cor- 
related interaction potentials of integer and noninteger indices. J. Math. Chem. 
49,290-295(2011) 

97. Guseinov, 1.1.: Unifed treatment of complex and real rotation-angular functions 
for two-center overlap integrals over arbitrary atomic orbitals. J. Math. Chem. 
49, 1011-1013(2011) 

98. Guseinov, I.I., Aksu, H.: Ground state energy calculations of isoelectronic series 
of He in double-zeta approximation using Coulomb potential with noninteger 
indices. Chin. Phys. Lett. 25, 896 - 898 (2008) 

99. Guseinov, 1. 1., Aydin, R., Bagci, A.: Application of complete orthonormal sets 
of f "-exponential-type orbitals to accurate ground and excited states calcu- 
lations of one-electron diatomic molecules using single-zeta approximation. 
Chin. Phys. Lett. 25, 2841 - 2844 (2008) 

100. Guseinov, 1.1., §ahin, F., Aydin, R., Bagci, A.: Use of basis sets of f "- 
exponential type orbitals in calculation of electronic energies for one-electron 
diatomic molecules by single-zeta approximation. Phys. Scr 77, 045,302-1 - 
045,302-6 (2008) 

101 . Guseinov, 1.1., Ertiirk, M.: Construction of different kinds of atomic and molec- 
ular orbitals using complete orthonormal sets of f "-ETO in single exponent 
approximation. Chin. Phys. Lett. 25, 2444 - 2447 (2008) 

102. Guseinov, 1.1., Ertiirk, M.: Application of combined Hartree-Fock-Roothaan 
theory to isoelectronic series of atoms using noninteger n-generalized expo- 
nential type orbitals. Commun. Math. Comput. Chem. (MATCH) 61, 603 - 613 
(2009) 

103. Guseinov, 1.1., Ertiirk, M.: Use of noninteger n-Slater type orbitals in combined 
Hartree-Fock-Roothaan theory for calculation of isoelectronic series of atoms 
Be to Ne. Int. J. Quantum Chem. 109, 176 - 184 (2009) 

104. Guseinov, 1.1., Ertiirk, M., §ahin, E., Aksu, H.: Calculations of isoelectronic 
series of He using noninteger n-Slater type orbitals in single and double zeta 
approximations. Chin. J. Chem. 26, 213 - 215 (2008) 

105. Guseinov, 1.1., Ertiirk, M., §ahin, E., Aksu, H., Bagci, A.: Calculation of nega- 
tive ions of B, C, N, O and F using noninteger n Slater type orbitals. J. Chin. 
Chem. Socc. 55, 303 - 306 (2008) 



60 

106. Guseinov, 1. 1., Erturk, M., Sahin, E.: Use of combined Hartree-Fock-Roothaan 
theory in evaluation of lowest states of K[Ar]4-s'^3d^ and Cr'^[Ar]4-s'^3d^ iso- 
electronic series over noninteger «-Slater type orbitals. Pramana J. Phys. 76, 
109-117(2011) 

107. Guseinov, I.I., Gorgiin, N.S.: Calculation of multicenter electric field gradient 
integrals over Slater-type orbitals using unsymmetrical one-range addition the- 
orems. J. Mol. Model. 17, 1517 - 1524 (2011) 

108. Guseinov, 1.1., Gorgun, N.S., Zaim, N.: Calculation of multicentre nuclear at- 
traction integrals over Slater-type orbitals using unsymmetrical one-range addi- 
tion theorems. Chin. Phys. B 19, 043,101-1 - 043,101-5 (2010) 

109. Guseinov, 1.1., Mamedov, B.A.: Computation of multicenter nuclear-attraction 
integrals of integer and noninteger n Slater orbitals using auxiliary functions. J. 
Theor. Comput. Chem. 1, 17 - 24 (2002) 

110. Guseinov, 1.1., Mamedov, B.A.: Evaluation of overlap integrals with integer and 
noninteger n Slater-type orbitals using auxiliary functions. J. Mol. Model. 8, 
272 - 276 (2002) 

111. Guseinov, 1.1., Mamedov, B. A.: Use of addition theorems in evaluation of multi- 
center nuclear-attraction and electron-repulsion integrals with integer and non- 
integer n Slater-type orbitals. Theor. Chem. Ace. 108, 21-26 (2002) 

1 12. Guseinov, 1.1., Mamedov, B.A.: Calculation of molecular electric and magnetic 
multipole moment integrals of integer and noninteger n Slater orbitals using 
overlap integrals. Int. J. Quantum Chem. 93, 9-19 (2003) 

113. Guseinov, 1.1., Mamedov, B.A.: Evaluation of multicenter electronic attraction, 
electric field and electric field gradient integrals with screened and nonscreened 
Coulomb potentials over integer and noninteger n Slater orbitals. J. Math. 
Chem. 36, 113-121 (2004) 

1 14. Guseinov, 1.1., Mamedov, B. A.: Evaluation of multicenter one-electron integrals 
of noninteger u screened Coulomb type potentials and their derivatives over 
noninteger n Slater orbitals. J. Chem. Phys. 121, 1649 - 1654 (2004) 

115. Guseinov, 1.1., Mamedov, B.A.: On evaluation of overlap integrals with nonin- 
teger principal quantum numbers. Comm. Theor Phys. 42, 753 - 756 (2004) 

1 16. Guseinov, 1.1., Mamedov, B.A.: Unified treatment of overlap integrals with inte- 
ger and noninteger n Slater-type orbitals using translational and rotational trans- 
formations for spherical harmonics. Can. J. Phys. 82, 205 - 21 1 (2004) 

117. Guseinov, 1.1., Mamedov, B. A.: Calculation of multicenter electronic attraction, 
electric field and electric field gradient integrals of Coulomb potential over in- 
teger and noninteger n Slater orbitals. J. Math. Chem. 37, 353 - 364 (2005) 

118. Guseinov, 1.1., Mamedov, B.A.: Evaluation of one- and two-electron multicenter 
integrals of Yukawa-like screened central and noncentral interaction potentials 
over Slater orbitals using addition theorems. Int. J. Mod. Phys. C 16, 837 - 842 
(2005) 

119. Guseinov, I.I., Mamedov, B.A.: Use of auxiliary functions in calculation of mul- 
ticenter electronic attraction, electric field and electric field gradient integrals of 
screened and nonscreened Coulomb potentials over noninteger n Slater orbitals. 
Chem. Phys. 312, 223 - 226 (2005) 



61 

120. Guseinov, I.I., Mamedov, B.A.: Use of recursion and analytical relations in eval- 
uation of hypergeometric functions arising in multicenter integrals with nonin- 
teger n Slater type orbitals. J. Math. Chem. 38, 51 1 - 517 (2005) 

121. Guseinov, 1.1., Mamedov, B.A.: Computation of three-center overlap integrals 
over noninteger n Slater type orbitals using f "-ETO. J. Theor. Comput. Chem. 
6,641-646(2007) 

122. Guseinov, 1.1., Mamedov, B.A.: Accurate evaluation of overlap integrals of 
Slater type orbitals with noninteger principal quantum numbers using complete 
orthonormal sets of f "-exponential type orbitals. J. Math. Chem. 43, 1527 - 
1532 (2008) 

123. Guseinov, 1.1., Mamedov, B.A.: Calculation of one-electron multicenter inte- 
grals of Slater type orbitals and Coulomb-Yukawa like correlated interaction 
potentials with integer and noninteger indices using unsymmetrical one-range 
addition theorems. Bull. Chem. Soc. Japan 83, 1047- 1051 (2010) 

124. Guseinov, 1.1., Mamedov, B. A.: Erratum to "the use of unsymmetrical one-range 
addition theorems of Slater type orbitals for the calculation of intermolecular 
Coulomb interaction energy" [Chem. Phys. Lett. 501 (201 1) 594]. Chem. Phys. 
Lett. 503, 185(2011) 

125. Guseinov, 1.1., Mamedov, B.A.: The use of unsymmetrical one-range addition 
theorems of Slater type orbitals for the calculation of intermolecular Coulomb 
interaction energy. Chem. Phys. Lett. 501, 594 - 597 (201 1) 

126. Guseinov, 1.1., Mamedov, B.A., Andic, Z., Cicek, S.: Use of unsymmetrical one- 
range addition theorems of Slater type orbitals in molecular electronic structure 
determination. J. Math. Chem. 45, 702 - 708 (2009) 

127. Guseinov, 1.1., Mamedov, B.A., Andig, Z.: Application of combined open shell 
Hartree-Fook-Roothaan theory to molecules using symmetrical one-range ad- 
dition theorems of Slater type orbitals. J. Math. Chem. 47, 295 - 304 (2010) 

128. Guseinov, 1.1., Mamedov, B.A., Orbay, M.: Calculation of three-center electric 
and magnetic multipole moment integrals using translation formulas for Slater- 
type orbitals. Theor. Chem. Ace. 104, 407 - 410 (2000) 

129. Guseinov, 1.1., Mamedov, B.A., Ozdogan, T, Orbay, M.: Calculation of mag- 
netic multipole moment integrals using translation formulas for Slater-type or- 
bitals. Pramana J. Phys. 53, 727 - 731 (1999) 

130. Guseinov, 1.1., Mamedov, B.A., Rzaeva, A.M.: Computation of molecular inte- 
grals over Slater-type orbitals. Vll. Calculation of multicenter molecular inte- 
grals by single-center expansion methods using different translation formulas. 
J. Mol. Struc. (Theochem) 544, 205 - 21 1 (2001) 

131. Guseinov, 1.1., Mamedov, B.A., Siinel, N.: Computation of molecular integrals 
over Slater-type orbitals. X. Calculation of overlap integrals with integer and 
noninteger n Slater orbitals using complete orthonormal sets of exponential 
functions. J. Mol. Struc. (Theochem) 593, 71-77 (2002) 

132. Guseinov, 1.1., Rzaeva, A.M., Mamedov, B.A., Orbay, M., Ozdogan, T, Oner, 
E: Computation of molecular integrals over Slater type orbitals. 11. Calculation 
of electric multipole moment integrals using translation formulas. J. Mol. Struc. 
(Theochem) 465, 7-9 (1999) 



62 

133. Guseinov, 1. 1., Sahin, E.: Evaluation of one-electron molecular integrals over 
complete orthonormal sets of f "-ETO using auxiliary functions. Int. J. Quan- 
tum Chem. 110, 1803 - 1808 (2010) 

134. Guseinov, 1.1., Sahin, E.: Evaluation of two-center Coulomb and hybrid inte- 
grals over complete orthonormal sets of f "-ETO using auxiliary functions. J. 
Mol. Model. 17, 851 - 856 (2011) 

135. Harris, F.E., Michels, H.H.: The evaluation of molecular integrals for Slater- 
type orbitals. Adv. Chem. Phys. 13, 205 - 266 (1967) 

136. Hartree, D.R.: The wave mechanics of an atom with a non-coulomb central 
field. Proc. Cambridge Phil. Soc. 24, 89 - 132 (1928) 

137. Higgins, J.R.: Completeness and Basis Properties of Sets of Special Functions. 
Cambridge U. P, Cambridge (1977) 

138. Hobson, E.W.: The Theory of Spherical and Ellipsoidal Harmonics. Chelsea, 
New York (1965). Originally published by Cambridge U. P., Cambridge (1931). 

139. Homeier, H.H.H., Weniger, E.J., Steinborn, E.O.: Simplified derivation of a 
one-range addition theorem of the Yukawa potential. Int. J. Quantum Chem. 
44,405-411(1992) 

140. Huzinaga, S.: Molecular integrals. Prog. Theor. Phys. Suppl. 40, 52 - 77 (1967) 

141. Kato, T.: On the eigenfunctions of many-particle systems in quantum mechan- 
ics. Commun. Pure Appl. Math. 10, 151 - 177 (1957) 

142. Kaufmann, K., Baumeister, W.: Single-centre expansion of Gaussian basis func- 
tions and the angular decomposition of their overlap integrals. J. Phys. B 22, 1 
-12(1989) 

143. Klahn, B.: Die Konvergenz des Ritz'schen Variationsverfahrens in der Quan- 
tenchemie. Ph.D. thesis, Mathematisch-Naturwissenschaftliche Fakultat der 
Georg-August-Universitat zu Gottingen, Gottingen (1975) 

144. Klahn, B.: Reviewof linear independence properties of infinite sets of functions 
used in quantum chemistry. Adv. Quantum Chem. 13, 155 - 209 (1981) 

145. Klahn, B., Bingel, W.A.: Completeness and linear independence of basis sets 
used in quantum chemistry. Int. J. Quantum Chem. 11, 943 - 957 (1977) 

146. Klahn, B., Bingel, W.A.: The convergence of the Rayleigh-Ritz method in quan- 
tum chemistry. 1. The criteria for convergence. Theor. Chim. Acta 44, 9-26 
(1977) 

147. Klahn, B., Bingel, W.A.: The convergence of the Rayleigh-Ritz method in quan- 
tum chemistry. 11. Investigation of the convergence for special systems of Slater, 
Gauss and two-electron functions. Theor. Chim. Acta 44, 27 - 43 (1977) 

148. Klahn, B., Morgan 111, J.D.: Rates of convergence of variational calculations 
and of expectation values. J. Chem. Phys. 81, 410 - 433 (1984) 

149. Kranz, H.H., Steinborn, E.O.: Implications and improvements of single-center 
expansions in molecules. Phys. Rev. A 25, 66 - 75 (1982) 

150. Le Guillou, J.C., Zinn-Justin, J. (eds.): Large-Order Behaviour of Perturbation 
Theory. North-Holland, Amsterdam (1990) 

151. Levin, D.: Development of non-linear transformations for improving conver- 
gence of sequences. Int. J. Comput. Math. B 3, 371 - 388 (1973) 

152. Li, J., Zang, W., Tian, J.: Simulation of Gaussian laser beams and electron dy- 
namics by Weniger transformation method. Opt. Expr. 17, 4959 - 4969 (2009) 



63 

153. Lowdin, P.O.: Quantum theory of cohesive properties of soHds. Adv. Phys. 5, 1 
-172(1956) 

154. Magnus, W., Oberhettinger, R, Soni, R.P: Formulas and Theorems for the Spe- 
cial Functions of Mathematical Physics. Springer- Verlag, New York (1966) 

155. Meyerson, M.D.: Every power series is a Taylor series. Amer. Math. Monthly 
88,51-52(1981) 

156. Niukkanen, A.W.: Fourier transforms of atomic orbitals. I. Reduction to four- 
dimensional harmonics and quadratic transformations. Int. J. Quantum Chem. 
25,941-955(1984) 

157. Novosadov, B.K.: Hydrogen-like atomic orbitals: Addition and expansion the- 
orems, integrals. Int. J. Quantum Chem. 24, 1 - 18. (1983) 

158. Nozaki, A.: How to detect divergence of some series with positive terms. Be- 
haviormetr. 15(23), 51-56 (1988) 

159. Olde Daalhuis, A.B.: Uniform asymptotic expansions for hypergeometric func- 
tions with large parameters I. Anal. Applic. 1, 1 1 1 - 120 (2003) 

160. Olde Daalhuis, A.B.: Uniform asymptotic expansions for hypergeometric func- 
tions with large parameters II. Anal. Applic. 1, 121 - 128 (2003) 

161. Olde Daalhuis, A.B.: Uniform asymptotic expansions for hypergeometric func- 
tions with large parameters III. Anal. Applic. pp. 199 - 210 (2010) 

162. Olver, F.W.J. : Asymptotics and Special Functions. A. K. Peters, Natick, Mass. 
(1997). Originally published by Academic Press, New York (1974). 

163. Olver, FW.J., Lozier, D.W., Boisvert, R.F, Clark, C.W. (eds.): NIST Handbook 
of Mathematical Functions. Cambridge U. P., Cambridge (2010) 

164. Prugovecki, E.: Quantum Mechanics in Hilbert Space, 2 edn. Academic Press, 
New York (1981) 

165. Rainville, E.D.: Special Functions. Chelsea, Bronx, New York (1971). Origi- 
nally published by Macmillan, New York (1960). 

166. Reed, M., Simon, B.: Methods of Modern Mathematical Physics I: Functional 
Analysis, 2 edn. Academic Press, New York (1980) 

167. Roothaan, C.C.J. : New developments in molecular orbital theory. Rev. Mod. 
Phys. 23,69-89(1951) 

168. Ruedenberg, K.: Bipolare Entwicklungen, Fouriertransformationen und 
Molekulare Mehrzentren-Integrale. Theor. Chim. Acta 7, 359 - 366 (1967) 

169. Sack, R.A.: Generalization of Laplace's expansion to arbitrary powers and func- 
tions of the distance between two points. J. Math. Phys. 5, 245 - 251 (1964) 

170. Safouhi, H.: Bessel, sine and cosine functions and extrapolation methods for 
computing molecular multi-center integrals. Numer. Algor. 54, 141 - 167 
(2010) 

171. Safouhi, H.: Integrals of the paramagnetic contribution in the relativistic calcu- 
lation of the shielding tensor. J. Math. Chem. 48, 601 - 616 (2010) 

172. Sansone, G.: Orthogonal Functions. Krieger, Huntington, NY (1977). Revised 
edition. Originally published by Wiley, New York (1959). 

173. Santos, F.D.: Finite range approximation in direct transfer reactions. Nucl. 
Phys. A 212, 341 - 364 (1973) 

174. Schwartz, L.: Theorie des Distributions, 2 edn. Hermann, Paris (1966) 



64 

175. Shavitt, I.: The Gaussian function in calculations of statistical mechanics and 
quantum mechanics. In: B. Alder, S. Fernbach, M. Rotenberg (eds.) Methods 
in Computational Physics Vol. 2. Quantum Mechanics, pp. 1 - 45. Academic 
Press, New York (1963) 

176. Shibuya, T.I., Wulfman, C.E.: Molecular orbitals in momentum space. Proc. 
Roy. Soc. A 286, 376 - 389 (1965) 

177. Silverstone, H.J.: Expansion about an arbitrary point of three-dimensional func- 
tions by the Fourier-transform convolution theorem. J. Chem. Phys. 47, 537 - 
540(1967) 

178. Slater, J.C: Atomic shielding constants. Phys. Rev. 36, 57 - 64 (1930) 

179. Slater, J.C: Analytic atomic wave functions. Phys. Rev. 42, 33 - 43 (1932) 

180. Smeyers, Y.G.: About evaluation of many-center molecular integrals. Theor. 
Chim. Acta 4, 452 - 459 (1966) 

181. Steinborn, E.O., Filter, E.: Translations of fields represented by spherical- 
harmonic expansions for molecular calculations. 11. Translations of powers of 
the length of the local vector Theor Chim. Acta 38, 261-271 (1975) 

182. Steinborn, E.O., Filter, E.: Translations of fields represented by spherical- 
harmonic expansions for molecular calculations. 111. Translations of reduced 
Bessel functions. Slater-type i-orbitals, and other functions. Theor. Chim. Acta 
38,273-281 (1975) 

183. Steinborn, E.O., Weniger, E.J.: Advantages of reduced Bessel functions as 
atomic orbitals: An application to Hj. Int. J. Quantum Chem. Symp. 11, 509 - 
516(1977) 

184. Steinborn, E.O., Weniger, E.J.: Sequence transformations for the efficient eval- 
uation of infinite series representations of some molecular integrals with expo- 
nentially decaying basis functions. J. Mol. Struct. (Theochem) 210, 71-78 
(1990) 

185. Stieltjes, T.J.: Recherches surquelques series semi-convergentes. Ann. Sci. Ec. 
Norm. Sup. 3, 201 - 258 (1886) 

186. Szego, G.: Orthogonal Polynomials. American Mathematical Society, Provi- 
dence, Rhode Island (1967) 

187. Temme, N.M.: Large parameter cases of the Gauss hypergeometric function. J. 
Comput. Appl. Math. 153, 441 - 462 (2003) 

188. Temme, N.M.: Numerical aspects of special functions. Acta Numer 16, 379 - 
478 (2007) 

189. Tricomi, E.G.: Vorlesungen iiber Orthogonalreihen, 2 edn. Springer- Verlag, 
Berlin (1970) 

190. Trivedi, H.P, Steinborn, E.O.: Numerical properties of a new translation for- 
mula for exponential-type functions and its application to one-electron multi- 
center integrals. Phys. Rev. A 25, 1 13 - 127 (1982) 

191. Weissbluth, M.: Atoms and Molecules. Academic Press, New York (1978) 

192. Weniger, E.J.: Untersuchung der Verwendbarkeit reduzierter Besselfunktionen 
als Basissatz fiir ab initio Rechnungen an Molekiilen. Vergleichende Rechnun- 
gen am Beispiel des Hj". Diplomarbeit, Fachbereich Chemie und Pharmazie, 
Universitat Regensburg (1977) 



65 

193. Weniger, EJ.: Reduzierte Bessel-Funktionen als LCAO-Basissatz: Analytische 
und numerische Untersuchungen. Ph.D. thesis, Fachbereich Chemie und Phar- 
mazie, Universitat Regensburg (1982). A short abstract of this thesis was pub- 
Hshed in Zentralblatt fur Mathematik 523, 444 (1984), abstract no. 65015 

94. Weniger, E.J.: Weakly convergent expansions of a plane wave and their use in 
Fourier integrals. J. Math. Phys. 26, 276 - 291 (1985) 

95. Weniger, E.J.: Nonlinear sequence transformations for the acceleration of con- 
vergence and the summation of divergent series. Comput. Phys. Rep. 10, 189 - 
371 (1989). Los Alamos Preprint math-ph/0306302 (http : //arXiv . org) 

96. Weniger, E.J.: On the summation of some divergent hypergeometric series and 
related perturbation expansions. J. Comput. Appl. Math. 32, 291 - 300 (1990) 

97. Weniger, E.J.: Interpolation between sequence transformations. Numer. Algor 
3,477-486(1992) 

98. Weniger, E.J.: Computation of the Whittaker function of the second kind by 
summing its divergent asymptotic series with the help of nonlinear sequence 
transformations. Comput. Phys. 10, 496 - 503 (1996) 

99. Weniger, E.J.: Addition theorems as three-dimensional Taylor expansions. Int. 
J. Quantum Chem. 76, 280 - 285 (2000) 

200. Weniger, E.J.: Prediction properties of Aitken's iterated A^ process, of Wynn's 
epsilon algorithm, and of Brezinski's iterated theta algorithm. J. Comput. Appl. 
Math. 122, 329 - 356 (2000). Reprinted in: C. Brezinski (ed.) Numerical Anal- 
ysis 2000, Vol. 2: Interpolation and Extrapolation, pp. 329 - 356. Elsevier, Am- 
sterdam (2000). 

201. Weniger, E.J.: Addition theorems as three-dimensional Taylor expansions. II. B 
functions and other exponentially decaying functions. Int. J. Quantum Chem. 
90,92-104(2002) 

202. Weniger, E.J.: The spherical tensor gradient operator. Collect. Czech. Chem. 
Commun. 70, 1225-1271 (2005) 

203. Weniger, E.J.: Asymptotic approximations to truncation errors of series repre- 
sentations for special functions. In: A. Iske, J. Levesley (eds.) Algorithms for 
Approximation, pp. 331 - 348. Springer- Verlag, Berlin (2007) 

204. Weniger, E.J.: Extended Comment on "One-Range Addition Theorems 
for Coulomb Interaction Potential and Its Derivatives" by I. I. Guseinov 
(Chem. Phys. Vol. 309 (2005), pp. 209 - 213). Los Alamos Preprint 
arXiv:0704.1088v3 [math-ph] (http: //arXiv. org) (2007) 

205. Weniger, E.J.: Further discussion of sequence transformation methods. 
Subtopic "Related Resources" (Rl) on the Numerical Recipes (Third Edition) 
Webnotes page http : //www . nr . com/webnotes/ (2007) 

206. Weniger, E.J.: Reply to "Extended Rejoinder to "Extended Comment on 
"One-Range Addition Theorems for Coulomb Interaction Potential and Its 
Derivatives" by I. I. Guseinov (Chem. Phys., Vol. 309 (2005), pp. 209-213)", 
arXiv:0706.0975v2". Los Alamos Preprint arXiv : 0707 . 3361vl [math-ph] 
(http : //arXiv . org) (2007) 

207. Weniger, E.J.: On the analyticity of Laguerre series. J. Phys. A 41, 425,207-1 
- 425,207^3 (2008) 



66 

208. Weniger, E.J.: The strange history of B functions or how theoretical chemists 
and mathematicians do (not) interact. Int. J. Quantum Chem. 109, 1706 - 1716 
(2009) 

209. Weniger, E.J.: An introduction to the topics presented at the conference "Ap- 
proximation and extrapolation of convergent and divergent sequences and se- 
ries" CIRM Luminy: September 28, 2009 - October 2, 2009. Appl. Numer. 
Math. 60, 1184-1187(2010) 

210. Weniger, E.J.: Summation of divergent power series by means of factorial se- 
ries. Appl. Numer. Math. 60, 1429-1441 (2010) 

211. Weniger, E.J.: One-range and two-range addition theorems. Topical review for 
Journal of Physics A, in preparation (201 1) 

212. Weniger, E.J., Cizek, J.: Rational approximations for the modified Bessel func- 
tion of the second kind. Comput. Phys. Commun. 59, 471 - 493 (1990) 

213. Weniger, E.J., Cizek, J., Vinette, P.: The summation of the ordinary and renor- 
malized perturbation series for the ground state energy of the quartic, sextic, 
and octic anharmonic oscillators using nonlinear sequence transformations. J. 
Math. Phys. 34, 571 - 609 (1993) 

214. Weniger, E.J., Kirtman, B.: Extrapolation methods for improving the con- 
vergence of oligomer calculations to the infinite chain limit of quasi- 
onedimensional stereoregular polymers. Comput. Math. Applic. 45, 189 - 215 
(2003) 

215. Weniger, E.J., Steinborn, E.O.: Programs for the coupling of spherical harmon- 
ics. Comput. Phys. Commun. 25, 149 - 157 (1982) 

216. Weniger, E.J., Steinborn, E.O.: The Fourier transforms of some exponential- 
type functions and their relevance to multicenter problems. J. Chem. Phys. 78, 
6121-6132(1983) 

217. Weniger, E.J., Steinborn, E.O.: Numerical properties of the convolution theo- 
rems of B functions. Phys. Rev. A 28, 2026 - 2041 (1983) 

218. Weniger, E.J., Steinborn, E.O.: Addition theorems for B functions and other 
exponentially decHning functions. J. Math. Phys. 30, 774 - 784 (1989) 

219. Wynn, P.: On a device for computing the e„t{Sn) transformation. Math. Tables 
Aids Comput. 10, 91 - 96 (1956) 

220. Zhang, S., Jin, J.: Computation of Special Functions. Wiley, New York (1996) 



